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ABSTRACT 

Our previous identification and spectroscopic confirmation of 431 faint, new planetary 
nebulae in the central 25 deg 2 region of the LMC permits us to now examine the shape 
of the LMC Planetary Nebula Luminosity Function (PNLF) through an unprecedented 
10 magnitude range. The majority of our newly discovered and previously known PNe 
were observed using the 2dF, multi-object fibre spectroscopy system on the 3.9-m 
Anglo- Australian Telescope and the FLAMES multi-object spectrograph on the 8-m 
VLT. We present reliable [Oiii]5007A and H/3 flux estimates based on calibrations to 
well established PN fluxes from previous surveys and spectroscopic standard stars. The 
bright cutoff (M*) of the PNLF is found by fitting a truncated exponential curve to 
the bright end of the PNLF over a 3.4 magnitude range. This cutoff is used to estimate 
a new distance modulus of 18.46 to the LMC, in close agreement with previous PNLF 
studies and the best estimates by other indicators. The bright end cutoff is robust to 
small samples of bright PNe since significantly increased PN samples do not change 
this fiducial. We then fit a truncated exponential curve directly to the bright end of the 
function over a 6 magnitude range and test the curve's ability to indicate the position 
of M*. Because of the significant increase in the number of LMC PN, the shape of 
the PNLF is now examined in greater detail than has previously been possible. Newly 
discovered features include a small increase in the number of PNe over the brightest 4 
magnitudes followed by a steep rise over 2 magnitudes, a peak at 6 magnitudes below 
the bright cutoff and an almost linear drop-off to the faint end. Dips at the bright 
end of the PNLF are examined in relation to the overall shape of the PNLF and 
the exponential increase in the number of PNe. Through cumulative functions, the 
new LMC PNLF is compared to those from the SMC and a new deep local Galactic 
sample revealing the effects of incompleteness. The new [Oiii]5007A LMC PNLF is 
then compared to our new H/3 LMC PNLF using calibrated and measured fluxes for 
the same objects, revealing the effects of metallicity on the [Oiii]5007A line. 

Key words: surveys - planetary nebulae: general - luminosity function - galaxies: 
Magellanic Clouds. 



1 INTRODUCTION 

The PN luminosity function (PNLF) describes the number 
density of PNe within a given system over any given luminos- 
ity range based on a prominent emission line such as [Om]. 
Its most frequent use is as an extragalactic distance indicator 
where it is the only standard candle that can be applied to 
all the large galaxies in the Local Supercluster (Jacoby, 1980, 
Ciardullo & Jacoby et al. 1989, Ciardullo & Jacoby, 1992, 
Ciardullo, 2006). In this role, the distance estimate relies on 
the shape created by the binning of [O ill] 5007A fluxes at the 
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bright end of the PNLF, which appears to be essentially in- 
variant across all galaxies tested thus far, regardless of their 
morphological type or metallicity (Mendez et al. 1993; Cia- 
rdullo et al. 2004) . There is a consistent [O ill] 5007A high 
luminosity cut-off beyond which no PNe are observed. The 
bright cut-off point is identified by fitting an exponential 
curve to the bright end of the PNLF and is believed to be in- 
variant for all PN systems in all galaxies (Stanghellini 1995, 
Ciardullo 2006). 

The concept of using PNe as a cosmological distance 
indicator was first suggested in the 1960's (Henize et al. 
1963; Hodge 1966). It was slow to gain acceptance because 
individual PNe are not standard candles, distances to local 
Galactic PNe were very difficult to estimate and extragalac- 
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tic PNe were severely under-sampled. The first study using 
an [O in] 5007A -emission-line based PN luminosity function 
as a distance indicator was performed in 1990 by Jacoby 
et al. (1990). Today, with large numbers of PNe being dis- 
covered in external galaxies (Kniazev et al. 2005, Ciardullo 
2006, Jacoby 2006) it is recognised as one of the most impor- 
tant and resilient distance indicators in extragalactic astron- 
omy (Ciardullo 2006). The PNLF method has been proven 
to be an accurate standard candle (Jacoby 1989; Ciardullo 
et al. 1989; Jacoby, Ciardullo, Ford, 1990; Jacoby, Walker, 
Ciardullo, 1990b) and appears to be extremely insensitive 
to parent stellar population (Ciardullo & Jacoby, 1992). 

The PNLF has other uses apart from distance determi- 
nation. Establishing the consistent shape of the function pro- 
vides a basis for estimating the number of PNe in any galaxy, 
given only the number of PNe in the brightest 3 or 4 magni- 
tudes. A PNLF which comprises a near complete sample of 
PNe across a galaxy can be used to estimate the luminosity- 
to-mass and dynamical age relations (Mendez et al. 1993). 
It has also been suggested as a valuable tool for studies of 
the initial-to-final mass relation and related mass loss pro- 
cesses in stars of early- type galaxies (Mendez et al. 1993). 
Lastly, the PNLF provides a unique probe into a galaxy's 
chemical and dynamical evolution. The mass and metallic- 
ity of the progenitor star largely determines the maximum 
[Om] line luminosity of the PN (Dopita et al. 1992). This 
makes the shape of the [O in] PNLF an important diagnostic 
for galactic chemical evolution. 

The first luminosity function for PNe in the LMC and 
SMC was constructed by Jacoby (1980). The advantage of 
using PNe in these local, external galaxies to create a lumi- 
nosity function was clear since the distance to all the PNe 
was essentially equal, faint PNe could still be identified and 
each galactic system could be studied in its entirety. The ma- 
jority of bright PNe in each galaxy had already been identi- 
fied so the bright end of the PNLF was able to be modeled 
and successfully used as a distance indicator. However, with 
a faint cutoff only 6 magnitudes below the brightest, the 
overall shape of the PNLF, especially at the faint end, was 
still unknown. 

The overall shape of the PNLF has previously been diffi- 
cult to determine for any galaxy due to sample incomplete- 
ness at the faint end. With 431 medium to faint PNe to- 
gether with 162 previously identified PNe now uncovered 
within the central 25 deg 2 of the LMC (Reid & Parker, 
2006b) [RPb] hereafter, we can now update the LMC PNLF 
and provide measured estimates to compare with previous 
theoretical simulations and other observed PNLFs. The RP 
PNe were discovered using an AAO/UKST Ha (+ effec- 
tively [Nil]) 40 field mini-survey of the entire LMC, SMC 
and surrounding regions. The details are provided in Reid & 
Parker (2006a) hence [RPa] and will not be repeated here. 
For details regarding the overall Ha survey, see Parker et 
al. (2005) and for the candidate selection technique and fol- 
lowup spectroscopic object confirmations, please see [RPb]. 

The new RP sample will enable both the shape and 
structure of the LMC PNLF to be analysed in detail for 
the first time over a much wider 10 magnitude range. The 
distance has always been estimated using an empirical (or 
theoretical) luminosity function and fitting this curve to the 
data using a statistically robust technique such as \ 2 or the 
maximum likelihood method to give the appropriate bright- 



end position. The bright-end intersection of this curve with 
the magnitude plane gives the position of M* (the brightest 
possible PN in the system) and the distance modulus. In 
this work, we use a cumulative function to determine the 
position of M* . We also fit the truncated exponential curve 
directly to the very complete bright end of the PNLF over 
a 6 magnitude range to test it's ability to depict PN evo- 
lution and indicate distance. Our ability to identify PNe to 
magnitudes as faint as 25 in [O ill] , 10 magnitudes below 
M*, gives us confidence that we are largely complete at the 
bright end of the PNLF. It also largely does away with the 
need to simulate and extrapolate the PNLF to account for 
unobserved, faint PNe with (M*-M) >5. 

A brief description of the PN spectra used in this study 
is provided in section [5] along with the data reduction pro- 
cedure. The new method of flux calibration for fibre spectra 
together with the de-reddening procedure is outlined in sec- 
tion |31 In section U we present the new [Om] PNLF for the 
LMC using 574 bright to faint PNe and compare the shape 
to the previous best LMC PNLF from Jacoby et al. (1990), 
and the empirical predictions of Ciadullo et al. (1989) and 
Mendez et al. (1993). In section [5] the PNLF is constructed 
using the H/3 line, a good predictor of central star tempera- 
ture, and compared to our new determination of the tradi- 
tional [Om]5007A based PNLF. 



2 2DF AND FLAMES OBSERVATIONS 
2.1 Brief description 

A five night observing run on the AAT using 2dF (Lewis et 
al. 2002) was undertaken in December 2004 to spectroscopi- 
cally confirm LMC emission candidates as PNe and to elim- 
inate contaminants such as Hn regions and emission line 
stars. 2dF was an ideal choice of instrument for the spec- 
troscopic followup of large numbers of candidate emission 
objects due to its unique ability to simultaneously observe 
400 targets (including objects, fiducial stars and sky posi- 
tions) with 2 arcsec fibres over a wide 2 degree diameter 
field area. The large corrector lens incorporates an atmo- 
spheric dispersion compensator, which is essential for wide 
wavelength coverage using small diameter fibres. 

The observations provided ~4,000 spectra. Individual 
exposure times were mostly 1200s using the 300B grat- 
ing with a central wavelength of 5852A and wavelength 
range 3600-8000A at a dispersion of 4.30A/pixel. These 
low-resolution observations, at 9.0A FWHM, were used as 
the primary means of object classification and provided the 
bright [O ill] 5007 fluxes for this study. All fields were then 
re-observed using the 1200R grating for kinematic studies 
[RPb] and improved resolution of the [S n] 6716,673lA lines 
for determination of electron densities. 

The 2dF raw data were processed using the AAO 2dF 
data reduction system, 2dFDR[]]. This software can probe the 
multifibre spectra creating the necessary calibration groups 
(eg. BIAS, DARK, FLAT, ARC etc). As a calibration expo- 
sure is reduced it is inserted into the appropriate group. 

The software was instructed to perform a subtraction 
of background scattered light prior to the extraction. The 

1 http:/ /www.aao.gov.au/AAO/2df/software.html#2dfdr 
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background is determined by fitting a function through the 
'dead' or unused fibres in the image. 

A bias frame was obtained for each observed field. The 
mean of the bias frame has the bias strip removed. This strip 
is subtracted from the data and the bias strip is trimmed 
from the data frame. A variance array is then derived from 
the data with values determined from photon statistics and 
readout noise. The resulting frame was then used throughout 
the reduction process. 

The data reduction system performs a wavelength cali- 
bration using the information from the spectrograph optical 
model. This is then refined using data from the arc lamp 
exposure. The lines in the arc lamp exposure are matched 
against a line list. A cubic fit for each fibre is then performed 
to the predicted and measured wavelengths of all lines which 
are non-blended, not too wide and not too weak. This fit is 
then used to further refine the wavelengths. 

To perform the sky subtraction, the data is first cor- 
rected for the relative fibre throughput, based on a through- 
put map derived from the dedicated sky fibres. The relative 
intensities of the skylines in the object data frame are used to 
determine the relative fibre throughput. This method saves 
time, as no offsetsky observations are required. The median 
sky spectrum was calculated from the median of all the sky 
spectra normalized by their mean flux. Each of these spectra 
are continuum subtracted, using a boxcar median smooth- 
ing with a 201 pixel box to define the continuum. A robust 
leastsquares fit of the counts was performed in the differen- 
tial or continuum subtracted data fibres versus the counts 
in the differential or continuum subtracted median sky. As- 
suming that the sky was the dominant source of emission, 
the slope of this fit gives the fibre throughput. The robust 
fit is especially important when dealing with faint objects. 
The sky fibre spectra in the data were then combined and 
subtracted from each fibre. Cosmic rays were rejected au- 
tomatically during the process of combining two or more 
observations. 

Our data includes spectroscopic observations of a sub- 
sample of 21 PNe in dense regions of the LMC main bar, 
undertaken using the multi-object fibre spectroscopic sys- 
tem, FLAMES (Pasquini 2002) on the VLT UT2 over three 
nights in December 2004. Gratings LR2 and LR3 allowed 
us to cover the most important optical diagnostic lines 
for PNe in the blue including [Oni]4363A, Hen4686A, 
H/3 and [O in] 4959 & 5007A. Grating LR6 covered the 
Ha, [Nil] 6548 and 6583A lines as well as the [Sn]6716 & 
673lA lines for electron densities. Details are provided in 
[RPa] and [RPb] but we will produce a separate paper ex- 
amining all 420 objects (including PNe, Hn regions, and 
WR stars) observed using FLAMES. 

The VLT FLAMES spectra were reduced using iraf 
tasks imred, specred and ccdred for multi-spec files. Cos- 
mic rays were rejected when combining frames. Using the 
weighted intensity of the continuum, the iraf scombine task 
was then used to combine the three different wavelength por- 
tions of the spectrum into one spectral image. 

2.2 Continuing Object followup 

Every PN in the LMC sample, both previously known and 
new, has been re-analysed in order to identify and exclude 
interlopers which may effect the PNLF. This is particularly 



important at the bright end for distance determination. Ha 
and off-band red images were examined alongside high and 
low resolution optical spectra. In addition we have searched 
the SAGE (Meixner et al. 2006), Spitzer maps and overlaid 
false colours at 3.6, 4.5, 5.8 and 8/im to assist verification. 
Finally, with mosaic radio maps of combined LMC data from 
ATCA and Parkes we were able to re-classify 3 ultra-bright 
'True' PNe as contaminants due to their strong radio fluxes 
>3mJy (see Filipovic et al. 2009). This is despite optical 
spectra which would otherwise strongly indicate a PN. These 
compact emission objects are now classified as Hn regions. 
For a full list of re-classified objects, please see Table [T] 

In order to provide independent checks of our 2dF data 
and create the most accurate LMC PN sample possible, we 
have undertaken a deep multi-wavelength study of PN can- 
didates selected from our original discovery list that exhib- 
ited unusual characteristics and/or were deemed suspicious. 
In particular, due to the possibility of minor positioning er- 
rors, we have made 78 longslit observations of previous 'pos- 
sible' and 'likely' PNe candidates in the RP sample using 
the MSSSO 2.3m telescope in order to firm-up their status. 
Using the double-beam spectrograph (DBS) on this tele- 
scope, the visible range (3200-9000A), is split by a dichroic 
at around 6000A and fed into two spectrographs, with red 
and blue optimized detectors. We used the 300B and 316R 
gratings to obtain a resolution of 4A for each arm of the 
DBS. Reduction and extraction of 2.3m spectra were per- 
formed using the standard iraf tasks imred, specred, cc- 
dred and Figaro's task bclean. One dimensional spectra 
were created and the background sky was subtracted. Final 
flux calibration used the standard stars, LTT7987, LTT9239, 
LTT2415 and LTT9491. Results are given in TableQ] where 
8 'likely' and 12 'possible' PNe have been re-classified. 



3 FLUX CALIBRATION AND 

DE-REDDENING OF THE 2DF FIBRE 
SPECTRA 

Decent flux calibration of the 2dF LMC PN spectra was re- 
quired in order to combine PN spectra from different 2dF 
fields and to make meaningful comparisons between fibre 
spectroscopy and long-slit observations of individual objects. 
A reliable flux calibration permits quick conversion to mag- 
nitudes for luminosity studies. It also permits a study into 
other nebula characteristics such as excitation class, tem- 
perature, density and mass. Since the majority of the spec- 
tra were gained by multi-fibre spectroscopy, it was neces- 
sary to verify the reliability of the measured line intensi- 
ties. A comparison of the same objects observed in two or 
more overlapping 2dF field plate exposures revealed vari- 
ations in line strengths of up to 35% from field to field, 
clearly rendering the standard reduced spectra unsuitable 
for combination in their raw state. A further comparison 
was made between combined 2dF line intensities and pub- 
lished fluxes from the Magellanic Cloud Planetary Nebulae 
spectral catalogue ([MCPN] thereafter) (Stanghellini et al. 
20020- The results, using all 65 MCPN which correspond 

1 Magellanic Cloud Planetary Nebula data is available from 
http://archive.stsci.edu/hst/mcpn/MCPN 
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Table 1. A list of 26 RP objects which have now been re-classified or had their status changed following new high resolution, higher S/N 
spectral observations and multi-wavelength analysis. ELS = Emission-line star; PN = Planetary Nebula; SNR = Supernova Remnant; 
Prob. refers to classification probability where T = True; L = Likely; P = Possible. Under the heading 'Reason', (1) refers to analysis 
using a 2.3m telescope optical spectrum, (2) refers to analysis using follow up observations on 2dF, 'IR' refers to false infrared colours 
at 3.6/im and 4.5/^m 5.8/xm and 8fim indicating Hn regions or hot stars, 'Radio' strong radio source (>3mjy). Please see the text 
(section 12. 2|l for more details. 



RP cat. 


d A ( h m s) 




Previous 


Previous 


New 


New 


Reason 


TV 


j zuuu 


Tonnn 
JzUUU 


i r^i 
1JJ 


Prob. 


ID 


Prob. 







n^ /in 9fi i s 
UO 4U zo.lo 


1C\ H7 1 

- / U 00 U / . 1 


PN' 
1 IN 


r r 
i 


H II region 


T 


wcaK [yj 1I1J DUU / {L) 


1U0 


n^ /in a x. or 
UO 4U 40. ZD 


vn oc r\R v 
- ( U zo UO. i 


PM 
1 IN 


r r 
i 


H II region 


1 


Large angular size, Radio, IR 


1 io 


UO O 1 If .oo 


7n 07 AQ 3 

-iu ui 4y.o 


PN 
1 IN 




oiN n 


T 


large size, strong [S II] lines (1) 


"3^9 

OOi 


nc; s xi i a 

UO Ol 14.00 


fif: na z.a fil 
-OO ZO 04.01 


PN 
1 IN 


p 




1 


I ni ,, ro TTT i Kf\r\7 (~\\ 

1 j o W 111J OUU 1 j np \*-) 


a on 
4yu 


n^ Q7 Q1 79 
UO O ( ol. / o 


71 in /IO 8 
- ( 1 1U 4o.o 


PINT 
1 IN 


p 


symbiotic 


i 


continuum peaking at 7500A(1) 


Oil 


n^ Q7 n« qq 
UO ol UO.oo 


A T 1 7 Q 
-Oy 4 ( 1 ( .O 


P M 
1 IN 


p 


H II region 


T 


[KJ 1I1J OUU 1 j rip — DU/o [1) 


fifi7 

OD t 


C\r. qn 9fi on 
UO ou zu.zu 


- I U IO U-L.O 


PN 


p 


ELS 


p 


Strong continuum (2) 


789 


98 18 sn 

UO ZO lO.OU 


-oy zo J-O.o 






PN 


L 


Good PN lines, morphology (2) 


798 


05 26 13.74 


-69 25 45.1 


PN 


L 


ELS 


L 


Strong continuum (1) 


841 


05 28 08.09 


-69 10 21.9 


PN 


P 


ELS 


L 


Strong continuum (1) 


856 


05 24 38.82 


-69 04 13.7 


PN 


T 


PN 


P 


[Om]5007/H/3 =160%, crowded environ. (2) 


872 


05 24 25.18 


-69 39 06.3 


PN 


L 


H II region 


T 


Strong continuum (1) 


993 


05 30 54.58 


-68 34 22.4 


PN 


L 


H II region 


T 


Strong radio source, IR 


1079 


05 30 33.09 


-66 57 41.5 


PN 


L 


ELS 


L 


Strong continuum (1) 


1113 


05 22 49.80 


-66 40 55.5 


PN 


P 


H II region 


T 


[O in] 5007 = H/3 (2) 


1192 


05 19 56.77 


-70 39 03.6 


PN 


T 


symbiotic 


L 


Red continuum (1) 


1434 


05 20 16.84 


-68 45 10.1 


PN 


L 


symbiotic 


L 


Red continuum (1) 


1495 


05 19 06.88 


-68 21 34.4 


PN 


P 


H II region 


T 


Radio 


1534 


05 21 29.65 


-67 51 06.6 


PN 


T 


PN + Hn 


L 


Radio, small PN 3" SE 


1541 


05 21 22.00 


-67 47 28.6 


PN 


L 


H II region 


P 


Radio, IR, [O m] 5007/IL3 = 60% 


1691 


05 00 32.23 


-70 00 49.0 


PN 


P 


ELS 


T 


Mostly H emission (2) 


1716 


04 54 24.49 


-69 29 42.5 


PN 


P 


H II region 


T 


Large size (2) 


1760 


05 03 51.47 


-68 57 23.7 


PN 


L 


ELS 


L 


Strong continuum (1) 


1783 


04 54 33.78 


-69 20 35.7 


PN 


P 


ELS 


T 


Low [Om]5007/H/3 (2) 


1933 


05 04 47.34 


-66 40 30.3 


PN 


P 


Hn region 


T 


Strong radio source (2) 


2194 


05 19 18.44 


-69 47 17.2 


PN 


P 


H II region 


T 


IR 



with the RP sample, are shown in Figure [T] General agree- 
ment is evident through an essentially log-linear relation but 
the scatter is too wide (sCl dex) to provide a reliable flux 
calibration. The variations in the [O in] 5007A line intensi- 
ties shown are mainly due to the different relative strengths 
in the [O ill] 5007A line from one 2dF field exposure to an- 
other in the raw, reddened spectra. A flux calibration of the 
spectra was required but to achieve this, a reliable method 
had to be found. This has proved quite difficult to do pre- 
viously with fibre-based spectra (eg. Mortlock et al. 2001; 
Georgakakis et al. 2004) . Below, an effective new method is 
described, which, in its careful application, is shown to be 
able to produce reliable flux calibration for the 2dF spectra. 

3.1 Flux calibration technique 

Altogether, 18 overlapping 2dF field-plate exposures were 
taken in order to cover the entire central 25deg 2 survey re- 
gion of the LMC. A comparison of line intensities from over- 
lapping areas revealed that each field plate/observational 
setup resulted in its own individual spectrum intensity scale. 
In addition, the 400 spectra from each field were fed to 2 
identical spectrographs. Although they were well matched, 
small instrumental differences were noticeable. Overall, the 
combined use of 2 different field-plates, each sending spec- 
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Figure 1. The uncorrected MCPN [O in] 5007A fluxes for 65 PNe 
are plotted against the [O III] 5007A 2dF 300B PN line intensities 
from all field exposures. General agreement is evident through an 
essentially log-linear relation with a turn-over at the bright end 
but the scatter is too wide (~1 dex) to provide a reliable flux 
calibration. The line represents a least squares fit to the data. 
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Figure 2. Four examples of the flux calibration procedure. In each case, an excellent agreement has been found between previously 
published [Om] and H/3 line fluxes and 2dF line intensities. The graphs on the left provide the calibration for CCDl and the graphs on 
the right provide the calibration for CCD2 from the same field-plate observation. Field CI is calibrated to 7 PNe, C2 to 6 PNe, HI to 
20 PNe and H2 to 28 PNe. The equation shown in the top left of each graph has been applied to all [O III] and H/3 line intensities for 
spectra fed to that particular CCD and field-plate combination. 



tra to 2 different spectrographs and CCD cameras resulted 
in the observed deviations seen when comparing previously 
published fluxes with line intensities for the same objects 
from multiple 2dF fields (Figure [1]). 

Experiments were undertaken in order to find the most 
reliable method of flux calibration using the observed line 
intensities. The best results were obtained by individually 
calibrating each spectral line on each field plate from each 
CCD camera to raw MCPN fluxes gained from HST ex- 
posures. The known PNe included and observed on each 
field plate were used as flux calibrators for each individ- 
ual measured field. In order to extend the reliability of 
the [Oiii]5007A flux calibration to magnitudes >19, the 
[Om]4959A and H/3486lA lines (Shaw et al. 2006) were in- 
cluded in the calibration. All three lines are close enough 
in the optical spectrum to produce reliable flux ratios. The 
individual H/3 and [O m]4959,5007A 2dF line intensities for 
known PNe observed for each ccd and each field plate ex- 
posure were graphed against the individual published fluxes 
for the same lines (eg. Figure [2]). The same method was ap- 



plied to our multi-fibre spectra from VLT FLAMES. All 2dF 
and FLAMES fields purposefully contained sufficient PNe 
with well calibrated fluxes to enable this process. The agree- 
ment of flux calibrated PNe from each spectrograph/field 
plate combination was considered strong enough (within 0.2 
dex, see Figure [2]) to allow its application to all the H/3 and 
[O m]4959,5007A emission lines for other PNe observed in 
the same field. In each line of best fit was derived 

and the underlying linear equation extracted. This equation 
became the calibrator for each emission line in each uncali- 
brated PN from the ccd and individual 2dF field plate expo- 
sure used. The equation was applied to each spectrum with 
a detectable line intensity in that field. 

Fluxes for LMC PNe from other catalogues (Jacoby et 
al. 1990, Leisy et al. private communication, Meathering- 
ham et al. 1988) were also included in order to build up the 
number of calibrators per field. Where only a de-reddened 
flux value was published, a relative reddening was applied 
in order to make a better correlation. Where 2 or more pub- 
lished fluxes were found for the same PN, stronger com- 
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parisons could be made. Each independent survey showed a 
mean agreement to within 0.2 dex, slightly increasing as the 
magnitude grew fainter. A comparison between calibrated 
[OlIl]5007A fluxes from 2dF and published fluxes from the 
MCPN, Jacoby et al. (1990) and Leisy et al. are shown in 
Figure [3] The agreement (<r< 0.2 dex) is very good across 
all surveys; especially since they represent space vs ground 
based observations in some cases with a variety of detectors 
and measurement methods. It also shows that the 2dF cali- 
brated fluxes remain within the a< 0.2 dex agreement. The 
mean slope of the fit (2dF/MCPN) is 1.01, where 1 would 
represent a perfect match. The median is exactly 1 and the 
standard deviation a = 0.12 dex. 

A spread in published line fluxes for the same PN 
is common but helps us to establish the uncertainties. 
Since the MCPN catalogue contained the largest number 
of PNe and the raw fluxes were 'as measured' (not de- 
reddened) these were given preference where irregularities 
became obvious. The MCPN set also includes some ground- 
based fluxes from ELCAT (Kaler et al. 1997) available at 
http: / / stsdas.s tsci.edu/elcat / 1 where the spectra are uncor- 
rected for extinction. The intention was to provide self con- 
sistent fluxes for all the 2dF observations. 

A calibrated quality and consistency check was then 
performed where fluxes for several PNe observed on more 
than one field plate exposure (due to the overlaps in 2dF 
pointings across the survey) were compared for consistency 
once the new flux calibration had been applied to each field. 
Figure [4] shows an excellent match in the high as well as the 
low flux regimes, proving that reliable and consistent fluxes 
can in fact be achieved this way. The equation to convert 
line intensities to fluxes in each field was then applied to the 
2dF and FLAMES [O ill] 5007A line intensities for each cor- 
responding field. Further line-to-flux calibration equations 
were then derived for other spectral lines such as Ha which, 
together with H/3, controls the de-reddening. A similar pro- 
cedure was then carried out for each of the spectral lines to 
provide their individual flux calibration. 

3.2 Reliability of the line intensities and fluxes 

At the distance of the LMC, most PNe are compact and 
point-like in ground-based observations allowing most of the 
flux to fall within the 2 arcsec diameter 2dF fibre. This is 
despite the larger reported angular diameters, which suffer 
from point spread function growth as a function of source 
intensity and include the very faint surrounding halos. The 
extended outer halos typically represent about 1,000th of the 
flux of the inner PN (Corradi, 2003). Clearly, the central re- 
gion of LMC PNe contains the majority of the light. There 
are 70 PNe out of our sample of 589 with inner shell diam- 
eters which exceed 2 arcsec by more than 1 arcsec. These 
are among the brightest LMC PNe and, despite the antic- 
ipated loss of light by 2dF spectroscopy, are still able to 
be accurately flux calibrated. For example, the loss of light 
from bright PNe is evidenced by the bright turnover in Fig- 
ure [T] However, individual spectrograph calibration which 
includes the brightest PNe in each field (Figure [2}, is able 
to largely correct for the assumed loss of light at the high- 
est magnitude level. It does this by calibrating the increasing 
brightness received by the fibre directly to a well established 
flux level. This level will vary for each spectrograph/plate 



observation. For this reason, calibration of 2dF spectra to 
a single standard star would produce erroneously low fluxes 
for bright and extended LMC objects. In the same way, if a 
longslit is too narrow, flux levels will be underestimated (Ja- 
coby & Kaler, 1993). This may be the cause of some small 
deviations in previously published fluxes for LMC PNe. 

Flux ratios of doublet or close lines from the central 
shell of the PN, such as the [S n] ratio used for the estima- 
tion of nebula densities, are largely consistent across the neb- 
ula and do not require flux calibration. The [S ll]6716,6731 
lines were measured from the high resolution 1200R spectra. 
All other lines, into the blue, were measured from the 300B 
spectra. 

3.2.1 Error estimations 

The line intensities measured from 2dF 300B and 1200R 
spectra have two main sources of error. The first is due to 
an uncertainty in where to place the continuum level for 
line intensity measurement. Since the author performed all 
the line measurements, they were all performed in the same 
manner, at the same base positions for any line, relative to 
the continuum and/or noise. Repeated measurements how- 
ever indicated an increased error estimate should apply to 
lines with a S/N of 5 or less. The error estimate increases 
from 10% at S/N = 5 to 40% at S/N = 3. About 5% of 
H/3 intensities are low enough that the increased measure- 
ment error should be applied. This is never applied to the 
[O III] 5007 line used for the PNLF in this study. These lines 
mostly have intensities >4 x H/3 and have repeat error esti- 
mates no greater than 3% for the faintest PNe. This system 
however provides a general error estimate for the line mea- 
surements in all 631 PNe. 

Secondly, there are difficulties in placing exact errors on 
the individual lines, due to the inherent variations in detec- 
tor response with wavelength. Luckily, this is the same for all 
objects on any given field plate exposure and can therefore 
be estimated from the response curve. The exposure times 
are also exactly the same across any given field and most of 
the exposures for the different fields are also of the same du- 
ration. With a small error included for these two variables, 
the remaining difficulty is the different absolute fluxes which 
cause a variation in the signal to noise. One fortunate as- 
pect of having a relatively large number of PNe re-observed 
on different 2dF fields is that a comparison between two or 
more observations of the same object is possible. A compar- 
ison of line intensity ratios for the same PN, observed on 
two or more fields (Figure [3]), gives an agreement to within 
±0.2 dex with ct=0.15 using a sample of 81 repeated PNe. 
This discrepancy includes ratios between strong lines and 
between strong to weak lines. 

In addition to the estimates for spectral line errors, the 
flux calibration adds an extra level of error. With a very 
good calibration of the 300B spectra to the MCPN fluxes 
(Figure [2}, the average difference to the published fluxes 
following calibration is ±0.18 dex with a — 0.5. These pub- 
lished fluxes will themselves have both measurement and 
systematic errors and these have been estimated at 0.2 dex 
through direct comparison. An allowance is made for these 
errors to increase with decreasing flux. This estimated in- 
crease is based on the standard deviation in the calibrated 
fluxes, increasing to ±0.3 a = 0.2 dex. 



No™ pm c i n +h e Large Magellanic Cloud (III): The PNLF 7 




-14.0 -13.5 -13.0 -12.5 -12.0 -11.5 -11.0 



Log 2dF & VLT calibrated [OIII]5007A fluxes (ergs cm 2 s 1 ) 

Figure 3. The calibrated [Om]5007A 2dF 300B line fluxes 
are plotted against the published [O III] 5007A fluxes from the 
MCPN on-line catalogue [MCPN], Jacoby, Walker. Ciadullo 
1990b, [JWC] and Leisy & Dennefeld 2006, [LD]. The sample of 
85 calibrated data points from 2dF and 13 from VLT FLAMES 
are drawn from 21 different fibre field exposures. Extremely good 
agreement can be seen between all the combined data points and 
published fluxes. Through this calibration method, 2dF and VLT 
FLAMES line intensities for all the remaining PNe observed in 
a given 2dF or FLAMES field exposure are reliably calibrated to 
fainter magnitudes. The mean agreement Aflux (RP-MCPN) is 
0.03±0.10 dex. A similar mean agreement is found between RP 
and JWC (-0.004±0.11) and RP and LD (0.02±0.13). 




2dF calibrated flux 

Figure 4. The [Om]5007A calibrated fluxes for 81 PNe with 
multiple observations (due to overlapping fields) are plotted in 
order to check the integrity of the calibration across different ob- 
servations including different CCD and field plate combinations. A 
good match (~0.2 dex) along the line of equality has been found, 
where there is a mixture of fields and flux intensities represented 
over a 4 dex range. 



With ground-based spectroscopy, atmospheric disper- 
sion can produce an extended object on the slit where blue 
light is separated from red light, also received by the detector 
(Filippenko, 1982). This effect can cause a systematic dif- 
ference when comparing standard stars observed at a lower 
air mass. It can also produce larger images in the blue than 



in the red due to more blue light being captured by the fi- 
bre. The increase is estimated as (1/A) ' 2 (Fried 1966; Boyd 
1978). To overcome this problem, 2dF has an 'Atmospheric 
Dispersion Compensator' (ADC) which largely compensates 
for this. With our flux calibration technique, however, lines 
from different regions of the optical spectrum are individ- 
ually calibrated to HST-based fluxes, largely circumventing 
problems of atmospheric dispersion and the differential air 
mass of standard stars. A correction for extinction is applied 
using the Balmer decrement, where the fixed ratio 1:2.86 ap- 
plies between Ha and H/3 in the absence of extinction. This 
corrected ratio is discussed in detail in the following subsec- 
tion. 

In addition to the photometric and flux calibration er- 
rors, an error estimate for the inclination of the LMC to the 
line of sight may be included. This could be as much as 0.04 
mag (Jacoby et al. 1990) based on a la depth uncertainty. 

One further problem can cause errors in flux estima- 
tions for LMC PNe. The LMC is filled with faint stars and 
diffuse emission of varying intensity. For individual slit spec- 
troscopy, it can be extremely difficult to find a region of 
blank sky along the slit. In the case of 2dF fibre spectroscopy, 
we dedicated ~40 fibres to specially selected sky positions. 
These positions represented average values of background 
sky in each field observed. The Ha map provided the most 
appropriate positions. For each field, the 2dfdr reduction 
program sampled and averaged these regions before apply- 
ing the sky subtraction to the objects. The resulting error es- 
timates for [O ill] 5007 as a function of line measurement and 
flux calibration are based on results from Figures [3] and [4] 
Line measurement errors increase from 2% to 10% and flux 
calibration errors increase from 6% to 20%. The line mea- 
surement error also includes an estimated error based on a 
function of the discrepancy between repeated, field to field 
and line ratio measurements. 

3.3 Corrections for extinction and reddening 

Extinction of light from distant objects is mainly the re- 
sult of interstellar dust. The light is both scattered and ab- 
sorbed, which increases the interstellar extinction towards 
shorter wavelengths. The amount of extinction will differ for 
each object and needs to be corrected to gain true fluxes. In 
the optical regime, the Ha/H/3 ratio was used to determine 
the extinction constant cH/3 (i.e., the logarithmic extinction 
at H/3 for each nebula. These hydrogen transitions are the 
strongest and easiest to accurately measure in the nebula 
spectrum and therefore provide a better estimate than other 
H lines. The observed Ha/H/3 ratio, when compared to the 
recombination value of 2.86 (Aller 1984), gives a logarithmic 
extinction at H/3 of: 

c(H/3) = (Zo ff (Ha/H/3) - %(2.86))/0.34 (1) 

This estimation is based on the relationship between 
observed and intrinsic intensities: 

I ob3 (H.a) = /mt(Ha) c( H/3)[/(Ha)-/(Hff)] (9 \ 
Iobs(m WH/3) ' [ > 

where[/ (Ha)-/ (H/3)] = -0.34 from the standard in- 
terstellar extinction curve given in Osterbrock (1989) and 
c(H/3) is the logarithmic extinction of H/3. The intrinsic ratio 
is mildly dependent on temperature: for T e between 2,500 
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and 20,000 K, the ratio varies from 3.30 to 2.76 (Osterbrock, 
2006). 

The value of c and 1(A) for all other lines A was then 
used to correct for interstellar extinction using the galactic 
extinction law from Whitford (1958) in the form of Miller & 
Mathews (1972): 



/cor (A) = / oi)s (A)10 



c[l+/(A)] 



(3) 



Reworking this equation for the H/3 line all other lines could 
be corrected for reddening using: 

Icor (A) _ / 6s(A) ^q C /(A) 



Icor(Hp) Iobs{Hp) 



(4) 



Application of the Balmer decrement ratio 2.86 between 
the Ha and H/3 lines is able to correct the spectrum in terms 
of the required ratios however the degree of c is dependent 
on flux calibration for each of these lines and any inter- 
nal inconsistencies. For fibre-based observations, the uncer- 
tainty in c(H/3) is extremely difficult to estimate as there 
may be small inconsistencies between fibres resulting from 
light transmission and position of the fibre on the plate. 
These will have an effect on the reddening corrected fluxes. 
Calibration of both the red and blue sides of the spectra each 
show a maximum uncertainty of 0.2 dex. This is consistent 
with the uncertainty found comparing published fluxes, in- 
dicating that internal inconsistencies are low. The reddening 
law of Howarth (1983) where E(B-V) = 0.689c was also em- 
ployed. It then follows that A5007 = 3.5 E(B-V) and A4861 
= 3.63 E(B-V). 

When estimating the errors, flux calibration is given 
predominant consideration together with the wide spectral 
range between the Ha and H/3 lines. This is estimated for 
Ha using the same method of 2dF flux calibration described 
above where the Ha spectral line in each nebula was individ- 
ually compared to published Ha fluxes for the same objects. 
The estimated uncertainty in c(H/3) as a result of combined 
maximum errors in line measurement and flux calibration 
estimations is <7%. 

The optical image of every PN was inspected for the 
presence of diffuse Hn across its immediate field. This was 
used as a means of testing the environment of PNe with high 
values of c(H/3) in the range 1.4 to 2.0. The available 8/im 
Spizer SAGE image (Meixner et al., 2006) was also exam- 
ined for areas where IR emission extended across individ- 
ual PNe. This was added confirmation that foreground dust 
could play a significant role in flux determination for some 
LMC PNe. The intensity of the ambient emission in Ha was 
then measured in counts, where, for a measure below 3 a, the 
noise begins to have a significant effect on the observed PN 
flux. Velocity measurements (RPb) for PNe and the velocity 
for the HI disk in the vicinity of each PN were then exam- 
ined to estimate whether the PN was foreground, within, 
or background to the bulk of observed emission. From the 
list of 292 newly discovered 'true' PNe, 11 were found in ar- 
eas of relatively strong H n emission. Of these, only 2 were 
identified as lying beyond the visible emission. Of the newly 
discovered PNe, all 4 'likely' PNe in areas of strong, dense 
Hn emission (>3 a the noise) exist within or beyond the 
Hi disk. All of these PNe have values of c(H/3) > 0.8. 

Due to the suspected low metallicity of many of the 
new, faint PNe in the RP LMC sample, as well as the gen- 
erally low metallicity of all PNe in the LMC compared to 
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Figure 5. The range of derived values for c(H/3) using the RP 
and previously known LMC PNe from 2dF 300B spectra. A peak 
at 0.4 shows the low overall extinction of most LMC PNe. Most 
of the bright, previously known PNe occupy the range 0-0.4. 



the Galaxy (Leisy & Dennefeld, 1996), we also consider the 
effect of collisional excitation on the Ha line (Stasinska, 
2002) and the effect it may have on the extinction calcu- 
lation. Because PNe are excited by power-law photoionisa- 
tion or shock-heating from the central star, collisional ex- 
citation and self-absorption can affect the intrinsic Balmer 
decrement. This can be caused by high T e resulting in the 
Ha/H/3 ratio no longer holding the canonical recombination 
value of 2.86 (Osterbrock & Ferland 2006). It is possible that 
this scenario applies to some PNe in the LMC. Since the de- 
reddening procedure involves Ha there may be some over- 
estimation of the extinction c(H/3). It has been shown that 
this error, which is independent of extinction, can be large 
when the lines have far separated wavelengths (Stasmska, 
2002). In low metallicity nebulae, there may be a range of 
evolutionary sequences where Ha/H/3 is significantly larger 
than the accepted recombination value. The resulting error 
has been estimated to be as high as 10%. Since there is no 
precise means to estimate the amount of correction required, 
care was taken to identify faint, low mass PNe with high T e 
and extremely high Ha/H/3 >10 ratios which may have been 
affected by power-law photoionisation or shock-heating (see 
Stasinska, 2002 fig. 2). In the LMC sample, 6 PNe with 
c(H/3) > 0.8 and T e > 25,000K have been identified. They 
are marked with a f in column 3 of appendix tables 3 and 4. 
Where inconsistencies (such as ultra- low N/O abundances) 
are found in 1 of these PNe (RP1584), the cause can be 
traced to collisional excitation or self absorption affecting 
c(H/3). 

The derived values for c(H/3) are shown in the histogram 
of Figure [5] Despite low extinction towards the LMC PNe 
we still find a few PNe with values of c reaching as high 
as 1.6. Most of the brighter previously known LMC PNe 
generally occupy the lower range from to 0.6 however some 
fainter Jacoby (1980) PNe in dusty, nebulous areas have 
c(H/3) values up to 1.0. Results for bright LMC PNe agree 
with values derived using fluxes in the MCPN catalogue. 
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The majority of these PNe have c(H/3) between and 0.6 
but also range as high as 1.6. For all the RP PNe, there is 
a peak at 0.4 with the largest proportion of PNe occupying 
the range 0.3 - 0.5. The mean extinction of 0.4 corresponds 
to E(B-V)— 0.31 mag, which is considerably higher than the 
line-of-sight reddening of 0.074-0.11 for stars in the LMC 
(Caldwell & Coulson 1986). It is also higher than a value 
of 0.2 found for LMC PNe by Herrmann & Ciardullo (2009) 
using a pre-RP catalogue sample to produce a basic skeleton 
of our Figure[5]plot. They cover less than half the magnitude 
range with lower sampling and base the Balmer decrement 
on estimated Ha values. Nonetheless, their sample of the 
bright end also shows evidence of an increase in extinction 
beginning at 3 magnitudes below the brightest. 

Somewhat higher values of c(H/3) in the RP sample are 
to be expected due to their faint magnitudes and often dusty 
environments. It does suggest that average internal dust red- 
dening may be as great as 0.5 mag. Most of the dust is ex- 
pected to consist of carbon-rich and dredged-up material, 
the quantity of which will have dependence on the mass and 
metallicity of the central star (Doptia et al. 1992). PNe at 
the bright end mainly occupy the medium excitation range 
(Reid & Parker 2006c) with central stars ranging in temper- 
ature from 80,000K to 126,000K making them very efficient 
at converting central star luminosity to luminosity in the 
nebula. Necessary corrections for c(H/3) were applied to the 
whole LMC sample on an individual object basis. The de- 
rived values of c(H/3) for each previously known and newly 
discovered PNe are given in appendix Tables 3 and 4 respec- 
tively. 

When deriving the distance to an extragalactic source 
using the PNLF, we use the raw [O in]5007A fluxes corrected 
only for foreground reddening. This is because none of the 
measured PNLFs and simulations to date have included the 
effects of circumstellar dust. Most surveys have been based 
on the [Om]5007A line. The H/3 line, required to derive c, 
is extremely difficult to measure from individual objects in 
external galaxies. Circumstellar dust therefore plays a large 
part in determining the observed brightness of the PNLF 
bright end. In order to make the correct comparison, our 
PNLF for distance determination has only been corrected 
for foreground reddening. The contribution of Galactic dust 
was removed assuming a foreground reddening value of E(B- 
V) = 0.074 (Caldwell & Coulson, 1985). 



4 THE NEW PNLF FOR THE LMC 

The PNLF is a plot of the total number of PNe in a defined 
volume at a particular magnitude at a given point of time. 
Since the [O ill] 5007A line is generally the brightest line op- 
tically emitted by the nebula, it has traditionally been used 
to plot the PNLF. The raw fluxes are listed in Appendix 
Tables 4 and 5 for the 164 previously known and 432 newly 
discovered PNe respectively. The [Oiii]5007A fluxes for the 
observed PN sample, corrected only for foreground redden- 
ing have been transformed to the magnitude system using 
the relation provided by Jacoby (1989): 

msooT = -2.5 log F5007 - 13.74. (5) 

Using the magnitudes, the luminosity function can be 
displayed through several methods. In principle we can get 
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Figure 8. In order to observe the effects of internal dust 
extinction within the LMC PNe, the foreground de-reddened 
[O m]5007A magnitudes are plotted against extinction. PNe in 
the brightest 4 magnitudes have low to medium extinction while 
PNe with medium to high extinction arc only found less than 4 
magnitudes below the brightest. These highly reddened faint PNe 
probably have intrinsically fainter cores which have evolved over 
a longer timescale. It is therefore possible that the dust properties 
of these objects may be rather different to those of the higher-core 
mass objects. 



an approximation to the luminosity function p(l) by plot- 
ting the number of PNe found within any number of equally 
spaced luminosity bins. Using this method, all previously 
known and new RP PNe within the central 25deg 2 of the 
LMC have been plotted into one of 54 bins, each 0.2 mag- 
nitude wide. This luminosity function, based on the fore- 
ground de-reddened flux of the [O in] 5007A emission line, 
in magnitudes, is shown as a histogram in Figure [6] The 
magnitude assigned to each bin represents the magnitude 
at the central position of each bin. Poisson error bars have 
been included. 

In the first PNLF derived for the Magellanic Clouds, Ja- 
coby (1980) using a sample of 41 PNe, demonstrated the the- 
oretical exponential function of Henize & Westerlund (1963), 
whereby a PN is treated as a uniformly expanding homoge- 
neous sphere ionised by a non-evolving central star. It fol- 
lows, therefore, that the number of PNe in each luminosity 
bin is proportional to the time those PNe spend at that 
luminosity (Jacoby, 1980, Ciardullo et al. 1989, Ciardullo 
et al. 2004). This may be represented by an exponential 
curve (Henize & Westerlund, 1963) with a sharp truncation 
designed to accommodate the bright end (Ciardullo et al. 
1989). The resulting curve is described by: 

Ar(M)oce°- 307M {l-e 3(M *- M) } (6) 

In this equation, the key parameter AT* is the absolute 
magnitude of the brightest possible PN. The limit to the 
high luminosity of the [O ill] 5007A line is mainly due to 
the initial mass of the progenitor star and its evolution to 
a white dwarf as a function of time in the ranges of 3-11 
Gyrs (Jacoby, 1989; Mendez et al., 1993; Stanghellini, 1995; 
Richer et al., 1997; Jacoby, 1997). This PNLF model has 
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Figure 6. The foreground de-reddened planetary nebulae luminosity function for the central 25 deg 2 of the LMC, with PNe separated 
into the categories of previously known and true, likely and possible RP PNe in the survey area for this study. The bins are 0.2 mags 
in width. Poisson error bars are shown. The errors are calculated individually for previously known, true, likely and possible PNe and 
combined for each bin. 
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Figure 7. The planetary nebulae luminosity function for the central 25 deg 2 of the LMC, with PNe separated into the categories 
of previously known and true, likely and possible RP PNe in the survey area for this study. All PNe are individually corrected for 
internal extinction using the Balmer decrement. The bins are 0.2 mags in width. Poisson error bars are shown. The errors are calculated 
individually for previously known, true, likely and possible PNe and combined for each bin. 
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been demonstrated as an excellent standard candle when 
compared to other distance indicators (Ciardullo, 2006). 

It has been shown that there is a weak correlation be- 
tween the luminosity of M* and the metallicity of a galaxy 
(Ciardullo & Jacoby, 1992; Dopita et al. 1992; Richer 1993; 
Ciardullo et al. 2002, 2005). This trend has been modeled 
by Dopita et al. (1992), showing that M* is bright for solar 
oxygen abundances but fades where metallicity increases or 
decreases away from solar. The correction for M* is given 
by: 

AM* = 0.928[O/H] 2 + 0.225[O/H] + 0.014 (7) 

where solar metallicity is assumed to be 12 + log (O/H) 
= 8.87 (Grevesse, Noels & Sauval 1996) (see Ciardullo et al. 
2002 for additional details). The metallicity of the LMC has 
been determined from the emission lines of Hll regions. A 
mean value of (12 + log [O/H]) 8.45 ±0.11 from LMC Hn re- 
gions was found by averaging mean values from Dufour, 
Shields & Talbot (1982), Dennefeld (1989), Monk, Barlow & 
Clegg (1988) and Ferrarese et al. (2000a). With this value, 
using the metallicity dependence of Doptia et al. (1992) we 
find a metal dependence (AM*) of 0.1 mag. 

In the LMC, the estimated metallicity is low compared 
to the Galaxy. The low oxygen abundance in LMC PNe, 
whether inherent or through depletion, also raises the elec- 
tron temperature thereby increasing the rate of collisional 
excitations per ion. Since oxygen is a coolant, a decrease in 
oxygen abundance only decreases the [O ill] 5007A flux by 
roughly the square root of the difference in abundance (Ja- 
coby 1989). The central star at the core of a PN affects the 
luminosity the opposite way to low metallicity. Low metal- 
licity in the progenitor star produces a more massive radius 
with greater UV flux (Lattanzio 1986). This energy compen- 
sates for the decreased emissivity of the nebula. The total 
[O ill] 5007A flux then becomes somewhat independent of 
metallicity. 

In order to reveal the effects of dust, internal to the 
LMC and the environment of each individual PN, we include 
another PNLF (see Figure [7]) where the same objects have 
been corrected for what we will refer to as 'internal extinc- 
tion' using the Balmer decrement. Both functions have been 
plotted on the same scale in order to clearly demonstrate 
the offset caused by extinction. A comparison of the fore- 
ground de-reddened (Figure [6]) and internally de-reddened 
(Figure [7| plots shows a highly consistent shape in the func- 
tion. 

Unfortunately, with most external galaxies, it is almost 
impossible to correct each observed PN for the effects of 
internal extinction. Remarkably however, Figure [7] implies 
that if the foreground de-reddened PNLF is simply shifted 
by applying 1 magnitude of extinction to each object, the in- 
ternally de-reddened PNLF is recovered, so that the shape of 
the observed and de-reddened PNLFs are essentially identi- 
cal. This result is rather surprising since although individual 
object extinctions vary between < cH/3 < 1.3 (0-2 mag), 
from a global perspective, the PNLF responds as though 
each PN has cH/3 = 0.4. The magnitude range and shape 
of the plot (Figure [7|, corrected for internal extinction, has 
only shifted ~1 magnitude brighter than that corrected for 
line of sight reddening. It implies that the one dimensional 
PNLF of external galaxies can be modeled by simply apply- 
ing a mean offset. 
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Figure 9. The planetary nebulae luminosity function for the cen- 
tral 25 deg 2 of the LMC, derived from an homogeneous sample of 
previously known and RP LMC PNe. In order to use the PNLF 
for distance determination, the data have only been corrected for 
line of sight reddening. The data have been binned into 0.2 mag 
intervals and plotted in log space. The solid line is the log of the 
truncated exponential curve (Equation 6) as predicted by Ciar- 
dullo et al. (1989) convolved for errors and placed at the best-fit 
position to the bright end of the observed function. Poisson error 
bars arc included. It is assumed that the decreasing number at 
magnitudes ni5007 > 21 is due to incompleteness. This provides 
the first, direct estimation of M* for [O III] 5007 in the LMC using 
a 3.4 magnitude range at the bright end of the LMC PNLF. 



In order to observe the effects of internal dust ex- 
tinction within the LMC PNe, the foreground de-reddened 
[Oiii]5007A magnitudes are plotted against extinction in 
Figure [8] PNe in the brightest 4 magnitudes have low to 
medium extinction (0 to 0.6) while PNe with medium to 
high extinction (>0.6) are only found > 4 magnitudes below 
the brightest. These highly reddened faint PNe probably 
have intrinsically fainter cores which have evolved over a 
longer timescale. It is therefore possible that the dust prop- 
erties of these objects may be rather different to those of the 
higher-core mass objects. Most of the PNe, across the full 
magnitude range, are located in the low to medium extinc- 
tion range. This suggests that, regardless of the diminishing 
intrinsic luminosity as a consequence of age, a large propor- 
tion of these PNe with high-mass cores have evolved fairly 
rapidly and their surrounding circumstellar envelope, which 
provides the extinction, have remained of similar opacity. 

This is not to say that there has not been evolution 
within each nebulae such as decreasing electron density 
against changes in mass. Although the trends shown in 
Figure [8] are strong, there are selection effects which may 
come into play. A further study including electron densities, 
masses and dynamical ages will be required in order to reveal 
these effects and understand the role played by circumstellar 
dust. 

4.1 New PNLF-based distance estimate to the 
LMC 

The standard means of determining the distance to a galaxy 
using the PNLF has been previously confined to the bright- 
est 1 or 2 magnitudes. We now extend this range to the 
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Table 2. Summary of the goodness of fit tests between the ob- 
served and empirical PNLFs. 



13 14 15 16 17 18 19 

Apparent [Olll] 5007A mag (corrected for internal extinction) 

Figure 10. The bright end of the planetary nebulae luminosity 
function to mag 19.3 for the central 25 deg 2 of the LMC, binned 
into 0.2 mag intervals and plotted in linear space. In order to 
examine the true shape of the [O m]5007-based PNLF, the data 
are fully corrected for internal extinction. The broken line is the 
theoretical truncated exponential curve (Equation 6) as predicted 
by Ciardullo et al. (1989). Error bars employ poisson statistics. 
The curve matches the bright end of the function very well. Data 
bins that exceed 2<r away from the theoretical curve are shown 
as open circles. This provides the first, direct estimation of the 
extinction corrected value of M* for [O III] 5007 in the LMC using 
the bright end of the LMC PNLF. 



brightest 3.4 magnitudes. With a good estimate of complete- 
ness to this magnitude depth, we avoid including the effects 
of any possible dip in the function, 4 magnitudes below the 
brightest. 

To estimate the distance, we use the x 2 method to fit 
a model PNLF to the observed PNLF corrected for fore- 
ground reddening (Figure [9|. The bright cut-off absolute 
magnitude (M* ) has been estimated as M* = -4.44 ± 0.05 
(Ciardullo et al. 1989, 2002; Jacoby et al. 1992). We adopt 
this value because it reflects the theoretical dependence of 
M* on metallicity (Dopita et al. 1992) and in order to make 
a direct comparison to previous LMC PNLF distance esti- 
mates using this value. This gives us a distance modulus 
of 18.46±0.2 for the LMC where the magnitude errors are 
used as systematic uncertainties for the distance modulus. 
The error estimate allows for errors in line measurement and 
calibration, a ler error in the fit, uncertainties in the distance 
calibrator galaxy, M31 and the shape model of the PNLF. 
The previous best PNLF for the LMC using real LMC PNe 
for a distance estimate is shown as Figure 4 in Jacoby et al 
(1990b). This figure includes the model exponential curve 
of Ciardullo et al. (1989). Figure [6] shows that we have only 
added one extra PNe to the bright end of the function. This 
being the case, it is not surprising that the distance modulus 
does not deviate greatly from the distance modulus of 18.44 
previously determined by Jacoby et al. (1990b). 

Speculation has arisen as to whether the shape of the 
PNLF itself can be used for distance determinations (eg. 
Mendez & Soffner, 1997) . Only a deep and complete survey 
can answer this question. 

Although it is extremely difficult to place an exact lim- 
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iting magnitude on our survey data, we are confident that 
we are complete to 7 mag below the brightest. We are, there- 
fore, in a good position to test if the truncated exponential 
curve (Ciardullo et al. 1989) is a good predictor of PNLF 
shape and can be used as a distance indicator. 

In order to test this, we fit the theoretical truncated 
exponential from Ciardullo et al. (1989), to the brightest 
6 magnitudes of the PNLF using the Levenberg-Marquardt 
fit method to achieve the best fit position. To check how 
well the observed PNLF follows the theoretical exponential 
curve (equation 6) we ran the Kolmogorov-Smirnov (KS), 
Anderson Darling and \ 2 tests since all three work very well 
with binned data. Three assumed distributions were tested 
in order to find the best fit. The results represent values of 
the 0.2 magnitude binning and no poor fitting data points 
have been smoothed or excluded. The KS test, based on 
the largest vertical difference between the theoretical and 
empirical cumulative distribution function, found a statistic 
of D=0.11 down the brightest 6 magnitudes assuming an 
exponential distribution. 

The threshold value of the significance level (P-value) of 
0.677 between the observed and theoretical functions means 
the null hypothesis (Ho) that the observed distribution does 
not match the theoretical distribution must be accepted with 
a 95% confidence limit. The major departure from the theo- 
retical curve occurs 1.5 magnitudes below the brightest PN 
and lasts for 1.5 magnitudes. This may be referred to as a 
dip in the observed function and will be discussed in subsec- 
tion l4.2l Despite peaks and troughs in the observed function, 
if one were to estimate PN populations simply by integrat- 
ing equation 6 over 6 magnitudes one would over estimate 
the population by only ~3%. 

The same objects, after correction for internal extinc- 
tion, are plotted in linear space (Figure I10[) in order to fit 
and test the truncated exponential curve while examining 
the bright end with increased detail. Error bars represent 
poisson statistics. The same truncated exponential curve, 
again has a mean bin-to-bin agreement of 87 ±53% against 
the observed data over 6 magnitudes. The goodness of fit 
tests between the observed and empirical PNLFs are sum- 
marised in Tabled Clearly, the exponential increase in the 
number of PNe to be found with decreasing magnitude, first 
suggested by Henize & Westerlund et al. (1963) is equally 
invalid once internal extinction is removed. The close statis- 
tical results are also further proof of the identical shape of 
the observed and internally corrected PNLFs, discussed in 
section 4. 

An up to date compendium of distance estimates to the 
LMC is shown in Table [3] These results represent a vari- 
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ety of calibration objects, methods and extinction estimates. 
The PNLF method has been shown to be a very robust 
method (Ciardullo et al. 2005) across all galaxies tested. 
Our confirmed distance estimate of 18.46 ±0.20 is in very 
good agreement with Cepheid and RR Lyrae distances de- 
rived from optical observations. It is also in close agreement 
with the previous PNLF distance estimate of 18.44 ±0.18 as 
found by Jacoby et al. (1990) where the underlying univer- 
sal luminosity function is used as a probability distribution 
function. The estimated distance was given by the abscissa 
of graphed maximum likelihood solutions. 

It is extremely encouraging that a direct fit to the bright 
end of the PNLF using our considerable new data over the 
first 6 magnitudes has produced a result in keeping with the 
previous LMC PNLF work using far fewer PNe. This impor- 
tant result shows that the bright end cut-off fit is robust to 
far less complete PN samples obtained for a galaxy and that 
its use as a standard candle is not compromised by modest 
sampling of the bright end of the PNLF. 

The 1990 IAU LMC distance consensus of 50 kpc was 
presented by Mould (1990). Using an absolute magnitude of 
-4.44 for M* in the LMC, the brightest PN would need to 
have an apparent magnitude of 14.05 to equal the 50 kpc 
distance. An increase in the value of M* would move the 
PNLF-based distance toward higher distance estimates. Our 
new PNLF distance indicator places the LMC at precisely 
49.2 kpc ±0.2 kpc, very close to 48.7 kpc ±2.32 kpc, which 
is the mean distance and mean error found from all the 
distance results shown in Table [3] 

4.2 The shape of the PNLF 

Since the PN population in the LMC is now one of the most 
complete, observed samples at fixed distance available, mod- 
eling of the PNLF can be undertaken using this population. 
Although this study has extended the luminosity range to 
far fainter limits than previously achievable, and is there- 
fore more complete, there still exists the possibility that the 
most luminous PNe are over-represented due to all surveys 
being flux-limited. The large number of PNe now available 
from the RP sample has a strong effect on the shape of the 
PNLF fainter than 5 absolute magnitudes below the bright- 
est. When examining the overall shape of the LMC PNLF, it 
is important to remember that the sample includes PNe that 
are both optically thick and optically thin, hydrogen burn- 
ers and helium burners, type I and type II PNe (Peimbert 
1978), young, middle-aged and old with different evolution- 
ary characteristics and morphologies (Vassiliadis & Wood, 
1994; Dopita et al. 1992). It is beyond the scope of this pa- 
per to undertake evolutionary modeling of the PNe at this 
stage, however, the shape of the PNLF can now be examined 
in detail. Unless specified, we use the PNLF, corrected for 
internal extinction when discussing the shape of the PNLF 
and any evolutionary models. 

The peak number of PNe occurs at magnitude 19, which 
is half a magnitude fainter than half the magnitude range 
of the sample. The distribution on both sides of the peak 
is rather different. The bright end is relatively flat for the 
initial 3.5 mag before a steep and almost constant rise to 
the peak. The faint end drops off far more gradually with 
only three detections in the faintest bin. 

The rapid rise at the bright end, from mag 17.5 (~4 mag 



below the brightest) is a new major feature. It has not been 
seen before in an observed sample of PNe in any galaxy but 
was previously predicted due to the mean age of 3-8 Gyr for 
stars in the LMC (Stryker & Butcher 1981; Frogel & Blanco 
1983; Hyland 1991). This places the largest number of PNe 
close to a luminosity of log (L/L©) ~ 3.65-3.75 in keeping 
with previous predictions (Dopita et al. 1992). 

This rise in the function coincides with the faint end of 
a sudden dip or inflection point found in some extragalactic 
studies (Ciardullo 2006). Importantly, we don't see a strong 
'dip' (also known as the 'Jacoby dip') hinted at in the earlier 
but far more incomplete LMC PNLF (see Figure [6] and [Tj) . 
Although there is a small decrease in the number of PNe 
between magnitudes 16.4 and 17.2 it is not a strong feature. 
Rather, the function is relatively flat until mag 17.2. The dip 
in the number of PNe occasionally seen at ~3.5-4 mag below 
the bright cut-off (eg. Jacoby & De Marco, 2002) has been 
attributed to the decline in luminosity of the central star as 
it starts descending the White Dwarf (WD) cooling track 
(eg. Frew 2008). Alternatively, it may indicate that the PN 
evolutionary phase can occur at different luminosity levels 
for stars of different progenitor masses where 4 magnitudes 
below the brightest is representative of the largest group. It 
may also be evidence for the bimodal luminosity function 
expected from post-AGB stars (Vassiliadis & Wood, 1994). 
Interestingly, the PNe populations sampled in star form- 
ing galaxies exhibit a broad dip at 1.5-2 magnitudes below 
the brightest. Both PNLFs (Figures [6] and [7] corrected for 
foreground and internal extinction respectively, show some 
convincing evidence for the existence of this dip. This is 
characteristic of galaxies and models in which the evolution 
of [O in] 5007A in PNe is governed by the rapid evolution of 
a high mass core. 

By graphing the PNLF into 0.4 magnitude bins as a 
smoothed frequency distribution (Figure fTTI) . it becomes eas- 
ier to identify any possible evolutionary dips in the shape 
of the PNLF. By including the truncated exponential curve 
we have a clearer impression of the dip between magnitudes 
15 and 16. Overall, however, we see a relatively flat distri- 
bution over the brightest 3 to 4 magnitudes followed by a 
steep rise in the function. It agrees very well with the theo- 
retically simulated PNLF created by Mendez et al. (2008), 
where there is a sudden rise at the bright end followed by 
a very gradual decline for 4 magnitudes and a steep rise 
over 2 magnitudes to a peak in the function. The reason for 
the function remaining relatively flat for 4 magnitudes may 
be related to the bright LMC PNe being a combination of 
H-burning and He-burning central stars with a small propor- 
tion of low mass H-burners. Figure 7 in Mendez at el. (2008) 
shows that the function steadily rises with no dip if the sim- 
ulation includes too many central stars with masses as low 
as O.55M0. It conversely dips too low if there are too many 
high mass central stars. The relatively steady distribution 
over the brightest 4 magnitudes of the LMC PNLF found 
through this study adds weight to probability that central 
star mass plays a major role in shaping the bright end; espe- 
cially if we accept that the steep rise ~4 mags in represents 
the point at which the central star starts its descent down 
the WD cooling track. 

The number of PNe identified in our deep Ha survey 
begins to decline at magnitudes msoo7 > 19. Is this due to 
incompleteness or is there a real turnover close to this point 
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Table 3. A 20 year compendium of distance estimates to the Large Magellanic Cloud including our new estimate using the LMC PNLF. 



Method 


Survey 


Distance 


Distance 






Modulus 


kpc 


B stars 


Shobbrock and Visvanathan (1987) 


18.3 ± 0.2 


45. 7 1 


Cepheids 


Feast and Walker (1987) 


18.47 ± 0.15 


49.4 


Cepheids 


Welch et al. (1987) 


18.57 ± 0.05 


51.7 


Cepheids 


Visvanathan (1989) 


18.42 ± 0.04 


48.3 


Cepheids 


Paturel et al. (1997) 


18.7 ± 0.02 


54.9 


Cepheids 


Luri and Torra et al. (1999) 


18.35 ± 0.13 


46.7 


Cepheids 


Bono et al. (2002) 


18.53 ± 0.08 


50.8 


Cepheids 


Keller and Wood (2006) 


18.54 ± 0.018 


51.0 


LT Eclipsing Binary Systems 


Pictrzynski et al. (2009) 


18.50 ± 3% 


50.1 


LPV & Miras 


Berereat, Knapik and Rutily (1998) 


18.50 ± 0.17 


50.1 


Miras 


Feast (1988) 


18.28 ± 0.6 


45. 3 1 


Model Atmospheres 


Eastman and Kirshner (1989) 


18.45 ± 0.28 


48.9 


MS fitting 


Chiosi and Pigatto (1986) 


18.5 ± 0.1 


50.1 


MS fitting 


Schommer et al. (1984) 


18.2 ± 0.2 


43.6 1 


MS fitting 


VandenBerg and Poll (1989) 


18.4 


47.8 


Novae 


Capaccioli et al. (1990) 


18.70 ± 0.2 


54.9 


O stars 


Conti, Garmany and Massey (1986) 


18.3 ± 0.3 


45. 7 1 


PNLF 


Jacoby, Walker and Ciardullo (1990) 


18.44 ± 0.18 


48.7 


PNLF 


Dopita, Jacoby, Vassiliadis (1992) 


18.37 ± 0.15 


47.2 


PNLF 


Reid and Parker (this work) (2009) 


18.46 ± 0.2 


49.2 


Red Clump stars 


Grocholski et al. (2007) 


18.40 ± 0.04 


47.8 


RR Lyraes 


Reid and Strugnell (1986) 


18.37 ± 0.15 


47.2 


RR Lyraes 


Walker and Mack (1987) 


18.44 ± 0.05 


48.7 


RR Lyraes 


Alcock, Alves, Axelrod et al. (2004) 


18.43 ± 1.6 


48.5 


RR Lyraes 


Catelan and Cortes (2008) 


18.44 ± 0.11 


48.7 
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Figure 11. The planetary nebulae luminosity function for the 
central 25 deg 2 of the LMC in 0.4 magnitude bins. The overall 
shape of the PNLF begins to resemble that simulated by Mcndcz 
et al (2008, Fig. 7). The simulation resulted in a smooth distribu- 
tion which may be close to the truth once all the bright PNe across 
the LMC are included in the PNLF. The truncated exponential 
curve (Ciadullo et al. 1989) has been included for comparison. Al- 
though we point out the broadly linear distribution at the faint 
end, no scientific conclusions can be drawn from this result as 
the faint end of the PNLF will suffer from incompleteness. It is 
however worth noting that our survey, to a depth of mag 25 in 
[O in] 5007, shows this linear decline. 



in the function? The decline begins rather steeply over the 
next 4 mag but clearly becomes increasingly incomplete in 
the faintest two magnitudes of the survey (mag 23 — > 25). 
If the theoretical and empirical PNLFs (Jacoby 1989; Cia- 
rdullo et al. 1989) were employed beyond the peak magni- 
tude of 19 as found in this work, they would predict the 
presence of a large number of PNe by magnitude 25. Are 
they there? Is there sufficient mass in the LMC to permit 
the co-existence of tens of thousands of faint PNe at this 
end of the function? 

It would be very surprising if we have missed detecting 
hundreds of PNe more than 5 magnitudes brighter than our 
detection limit. Our PNe were selected using a deep Ha fil- 
ter (see RPa) which includes the [Nll]6548 and 6583A lines. 
Deep PN surveys of the LMC have already been conducted 
using [O in] filters with enormous success (RPa) . Our survey 
using Ha and [N n] has allowed us to discover a large num- 
ber of PNe with extremely low [O ill] levels. Where [O ill] is 
low in abundance, [Nil] is normally enhanced (RPa,b), per- 
mitting us to uncover a large population of PNe that would 
be very difficult to discover in an [O ill] survey. Even so, it 
is always possible that a deeper [O ill] survey and new ex- 
traction techniques may uncover more PNe. The only really 
hard constraint on the number of PNe to be found in the 
LMC comes from stellar evolutionary theory. The number 
of PNe shouldn't be greater than the product of the stellar 
evolutionary flux (2E-11 stars yr _1 L© -1 ) times the LMC 
bolometric luminosity times the PN lifetime. Our previous 
work (RPb) estimates the existence of 956 ± 141 PNe in the 
LMC, in close agreement with an earlier estimation of 996 
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by Jacoby (1980). The derived luminosity-specific number 
of PNe (pi) in the LMC according to Buzzoni et al. (2006), 
is given as log a = -6.57 ± 0.04, representing 1040 ± 60 PNe. 

We therefore suggest that a turnover is likely to occur 
at some midpoint in the function. At this stage, we can 
safely say that equation 6 represents the bright end of the 
PNLF with a range of 6 magnitudes. After this initial 6 
magnitudes and the peak in the PNLF, there is likely to be 
a flattening or even a decline in the number of PNe to be 
found. Only further deep surveys will be able to solve this 
issue conclusively. 



4.3 Luminosity implied time 

If we accept that the number of PNe in each luminosity bin 
is generally proportional to the time those PNe spend at 
that luminosity (Jacoby 1980; Ciardullo et al. 1989; Ciar- 
dullo et al. 2004), then we may conclude that PNe evolve 
away from the brightest 3 magnitudes at a proportionally 
fast pace. The time spent at each luminosity bin is plot- 
ted in Figure [T2"l where each luminosity bin is 1 magnitude 
wide. This histogram ignores the possible incompleteness of 
the faint end. It is also insensitive to any intrinsic varia- 
tion and 'birth-to-death' luminosity range for a given PN 
as progenitors range between 1 and 8 M©. Many PNe may 
have low-luminosity cores, and therefore join the luminosity 
function at fainter magnitudes. A simulated increase of the 
faintest 2 magnitudes would further decrease the percent- 
ages at the bright end. All this plot can effectively show is 
the number of PNe that may be found per luminosity bin 
at a single instant of time. 

Each bin has been scaled as a percentage of the number 
of PNe in the whole sample. Figure \TE\ shows that the per- 
centage of PNe in the brightest 4 magnitudes accounts for 
only 14.6% of all LMC PNe within the deep sample. In other 
words, 40% of the magnitude range contains only 14.6% of 
the PNe. It has been suggested that this drop in number 
density brighter than mag 18 is due to the negative slope of 
the initial mass function, especially effecting young stellar 
populations (Shaw, 1989) and the rapidly accelerating rate 
of evolution across the H-R diagram with increasing lumi- 
nosity, particularly evident on the hydrogen-burning tracks 
(Dopita et al. 1992). Simulations conducted by Mendez et 
al. (1993) suggest that these PNe at the bright end are also 
predominantly optically thin. 



4.4 The cumulative PNLF and comparisons with 
other galaxies 

In order to compare the depth of our derived PNLF with pre- 
vious LMC PNLF models and PNLF in other galaxies, we 
create a cumulative distribution of the function. We present 
the cumulative function as a percentage of PNe found at 
each magnitude and compare this to model cumulative func- 
tions. This approach helps to reveal the positions at which 
incompleteness begins to set in. It also shows the way in 
which incompleteness affects the overall shape of the cumu- 
lative function. 

The cumulative LMC PNLF for [O in] 5007 is shown in 
Figure [13] where it is compared with the cumulative PNLFs 
calculated in the same way for the SMC, using data from 
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Figure 12. The percentage of PNe at each luminosity bin rep- 
resents the percentage of time PNe spend at that luminosity ac- 
cording to predictions of Jacoby (1980); Ciardullo et al. (1989) 
and Ciardullo et al. (2004). This plot ignores the possible incom- 
pleteness of the faint end. It is also insensitive to any intrinsic 
variation and 'birth-to-death' luminosity range for a given PN as 
progenitors range between 1 and 8 Mq. A simulated increase of 
the faintest 2 magnitudes would further decrease the percentages 
at the bright end. 



Jacoby (2006) and for local PNe using the new highly com- 
plete local Galactic volume catalogue of Frew (2008) which 
includes PNe out to 3kpc. These objects are included in 
the MASH catalogue of Parker et al (2006) . Although each 
PNLF is affected by incompleteness, compared to the LMC, 
the local and SMC samples have a larger percentage of PNe 
in the brightest 6 magnitude range (41% for the LMC, com- 
pared to 54% for the SMC and 55% for the local PNLFs). 
The shape of the bright end of the PNLF is therefore very 
sensitive to the relative number of PNe included at the very 
brightest magnitude. It is also sensitive to the peak num- 
ber density, the flux limit of the survey and the decreasing 
number of PNe at the faint end. 

The large increase in the number of LMC PNe between 
4 and the peak at 7 magnitudes below the brightest causes 
the LMC PNLF to rise away from the predicted exponential 
curve. With an identical magnitude range, the local PNLF 
shows a more constant rise. The SMC PNLF, with its smaller 
magnitude range, shows a small dip between 4 and 5 mag 
below the brightest, after which it also rises steeply. The 
steepest section of each cumulative function represents a rise 
to the peak in number density. This feature is easiest to see 
in the LMC PNLF, however, it's presence also in the SMC 
and between 6 and 7.5 magnitude ranges in the local PNLF 
suggests that it may be common to PN evolution in most 
galaxies. 

The final three magnitudes of the LMC PNLF show 
a gradual turnover towards the faint drop-off. A stronger 
turnover is also seen in the local sample, affecting the fi- 
nal 2 magnitudes, with only 8% of PNe found in this range, 
prior to the faint drop-off. These turnovers are the result 
of incompleteness at the faint end, where we would expect 
the number of PNe to either rise or remain stable beyond 
the peak density (~6 mag below the brightest). We do not 
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discount the possibility that a certain drop in number den- 
sity towards the faint end may be a real feature suggesting 
that most PNe spend the majority of their lifecycle between 
4 and 7 magnitudes below the brightest. The SMC PNLF 
only shows a turnover (drop in number density) in the final 
magnitude of its faint end. This short turnover would nor- 
mally indicate that the SMC PNLF is more complete up to 
its final magnitude drop-off than the other two PNLFs. The 
SMC faint turnover, however, also coincides with the drop- 
off from the peak number density in the LMC PNLF. The 
LMC contains a further 228 PNe fainter than this magni- 
tude. Therefore, the faint end of the SMC PNLF is actually 
less complete than the other two and we predict more PNe 
be found >4 mag below the brightest and 3 magnitudes 
deeper than currently available. 

Using equation 6, Mendez et al. (1993) have obtained 
the cumulative PNLF for any magnitude M fainter than the 
cutoff by integrating N(M) between M* and M using the 
following formula: 

A „ 'AM* 

K(M) = Ife^-e^*]-^ [ e Mc 2 -3)_ e M*(c 2 - 3)](8) 

c z (c 2 — 3) 

The dotted lines in Figure [13] represent this predicted 
cumulative PNLF from Eq. (8), adjusted to fit the magni- 
tude range of PNe in the SMC and LMC. Comparing this 
predicted curve to the SMC PNLF suggests there is a lack 
of PNe in the range 4-6 below the brightest. This is affect- 
ing the overall shape of this PNLF and creating the upward 
'hump' from 1 to 4 mag below the brightest. The cumulative 
curve from Eq. (8), when fitted to the LMC, agrees very well 
at the brightest 4 magnitudes. After this, it does not follow 
the strong peak in the function but continues at a steady 
exponential rate. A larger number of PNe at the faintest 3 
magnitudes would effectively bring the LMC PNLF closer to 
the predicted curve. The local PNLF also shares a 10 mag 
range in [O ill] and so shares the same predicted curve as 
our LMC PNLF, but traces out a path much further to the 
bright end and away from the simulated PNLF. This posi- 
tion for the local PNLF, when compared to the LMC PNLF 
either suggests that there should be a great many more faint 
PNe (>6 mag below the brightest) to be found within the 
local 3 kpc radii range or the local PNLF is affected by selec- 
tion effects since it is only sampling a relatively small region 
of our Galaxy. 

While the formula of Ciardullo et al. (1989) results in a 
constantly increasing PNLF towards fainter magnitudes, the 
simulated PNLF of Mendez et al. (1997) produces two rises 
and declines. The first decline is between [O ill] 5007A mag- 
nitudes -3.5 and -2.5 (Mendez et al. 1997, fig. 4). It then 
rises again from mags fainter than -2.3 until it reaches a 
peak at ~-0.2. From this point, it begins a steady decline. 
This simulation however was re-worked to allow for PNe on 
cooling tracks by including a random distribution of the ab- 
sorbing factor between 0.1 and 1. With this adjustment, the 
simulation fits the formula of Ciardullo et al. (1989) very 
well (see Mendez et al. 1997, fig 5). This was a prediction, 
where severe incompleteness hampered any true test with 
observations. The simulation of Mendez et al. (1997, fig. 5) 
after a correction to allow absorbing factors to tend to 0, 
shows a peak in the function at mag ~6 and a faint end 
decline from that point. This is exactly what is seen in the 
LMC PNLF using the RP sample (Figure [T3l) . This in turn 




[OIII]5007A magnitudes below the brightest 

Figure 13. The cumulative PNLFs for [Om]5007 using the new 
PNLF for the central LMC (solid line - this work), the revised 
deep SMC PNLF of Jacoby (2006) and a deep local sample from 
Frew (2008), to a volume radii of 3.0kpc, centred on the sun. The 
dotted lines are the plots of Eq. (8) from Mendez et al. (1993) 
with ci = 12 for the SMC, ci = 6 for the LMC and local PNLF 
to 3kpc, C2 = 0.307 for both. 



strengthens the case for a strong rise in the function from 4 
magnitudes below the brightest. This rise is the result of an 
evolutionary change where a young, high mass and luminous 
PN reaches a peak in its photoionisation potential. It then 
shifts towards the WD cooling track with increasingly lower 
density to mass ratios. The sudden rise in number density 
shows that this change occurs rapidly. 

The extension of the PNLF to fainter magnitudes de- 
scribes complete population effects as predicted by the 
Shklovski model (1981) and permits us to observe what ef- 
fect the faint PN population may reveal for successful PNLF 
modeling. It is hoped that this will in turn help to further 
refine PNLF simulations. 



5 THE POTENTIAL OF THE LMC H/3 PNLF 

The PNLF is essentially based on the evolution of the central 
ionising star (Jacoby, 1989). In particular, the H/3 intensity 
is a measure of the central star's ionising luminosity (Oster- 
brock, 1989) and can also indicate the heating of the nebula 
as long as it is optically thick. This would indicate that a 
hydrogen emission line would also be a natural choice for a 
PNLF. Since the evolution of the central star's luminosity 
gives rise to the PN luminosity function, there should be a 
correlation between the relative strengths of the integrated 
[O ill] 5007A and H/3 fluxes. The H/3 line is also less influ- 
enced by the metallicity of the parent galaxy, which could 
make it an ideal choice for distance determination. A PNLF 
was constructed using H/3 in order to discover the bright 
cutoff and shape of the function. A comparison with the 
[O ill] 5007 PNLF should then reveal broad influences of the 
LMC's low metallicity. The H/3 flux by comparison is typi- 
cally > 5 times fainter than [O in] 5007 in most PNe making 
it much more difficult to spectroscopically measure in ex- 
ternal galaxies. This is why H/3 is not generally used for 
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Figure 14. The H/3 planetary nebulae luminosity function for the central 25 deg 2 of the LMC, with PNe separated into the categories of 
previously known, true, likely and possible PNe in the survey area for this study. Only foreground de-reddening has been applied to the 
raw magnitudes. The bins are 0.2 mags in width. Poisson error bars are included. The errors are calculated individually for previously 
known, true, likely and possible PNe as shown in Tables l4l &: |51 and then combined for each magnitude bin. 
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Figure 15. The H/3 planetary nebulae luminosity function for the central 25 deg 2 of the LMC, with PNe separated into the categories of 
previously known, true, likely and possible PNe in the survey area for this study. A correction for internal extinction has been applied to 
each PNe using the Balmer decrement. The bins are 0.2 mags in width. Poisson error bars have been included. The errors are calculated 
individually for previously known, true, likely and possible PNe as shown in Tables f4] & \E\ and then combined for each magnitude bin. 
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distance determination. H/3 has a bright cut-off at about 2.4 
mag fainter than that of [0 in] 5007A. 

For optically thick PNe, the H/3 flux should be a reliable 
indicator of stellar luminosity. The maximum conversion ef- 
ficiency from stellar luminosity to luminosity in the H/3 line 
only varies by a factor of 2, even though a maximum of just 
under 1% of stellar luminosity is converted into H/3 photons 
(Dopita et al. 1992). Dopita et al. (1992) have found the 
variation in this conversion efficiency to be even less. 

In order to convert our H/3 fluxes, expressed in ergs 
cm~ 2 s" 1 , to magnitudes, we use the definition given by 
Jacoby (1989): 

m f = -2.5 log / - 13.74 (9) 

where / represents the flux to be roughly converted to 
the broad V-band magnitude. The close interval between 
the [O iii]5007A and H/3 lines allows use of the same defini- 
tion to convert fluxes from ergs cm -2 s _1 , to magnitudes for 
both lines. This technique, previously adopted by Mendez 
et al. (1993), allows us to make a direct comparison between 
[Om]5007A and H/3 magnitudes on the same scale. The lu- 
minosity function for the H/3 line is shown in Figure [T4l for 
foreground de-reddened objects where the brightest PN is 
found at an apparent magnitude of m486i = 16.8 In Fig- 
ure [15] the PNLF is shown following correction for internal 
extinction where the brightest PN is found at an apparent 
magnitude of 16.2. 

In Figure [16] we show the H/3 PNLF as a log plot with 
the truncated exponential curve of Ciardullo et al. (1989) 
added for comparison. The curve has been placed at the 
best fit to the bright end, using the Levenberg-Marquardt 
fit method. Using the bright intersection of the curve with 
the s-axis, a value of M4 861 = 16.8±0.21 was found. This cor- 
responds to an absolute magnitude of -1.68±0.21 for MJggi. 
Statistically, the goodness of fit tests to the theoretical ex- 
ponential curve give F=1.31, a KS statistic of 0.13 and \ 2 
statistic of 2.85 providing 95% confidence limits that the 
observed PNLF is unlikely to be drawn from the empirical 
one. At present there is little data on extragalactic H/3 fluxes 
available with which to compare these results. Stanghellini 
(1995) constructed an Rf3 PNLF for the Magellanic Clouds 
based on the fluxes of Richer (private communication) to- 
gether with a simulation of 1000 optically thick PNe using 
Salpeter's IMF and Weidemann's initial mass-final mass re- 
lation, with M max = 0.7 M (Weidemann 1987). This data 
produced very encouraging comparisons, however, variations 
in the mass truncation produced significantly different re- 
sults. Using M m ax = 0.7 Mg, the peak in the distribution 
occurred only 1 magnitude below the brightest PN. Increas- 
ing Mmax to 1.4 Mq shifted the peak 2.2 magnitudes fainter 
but the entire distribution only covered a 4 magnitude range. 

The bright end only is plotted in linear space in Fig- 
ure [17] using the same data corrected for internal extinction. 
Analogous to the [Oin]5007A function, the brightest PN in 
the H/3 function also increases 0.8 magnitudes after correc- 
tion. Assuming an exponential growth in the function, and 
comparing the empirical PNLF of Ciardullo et al. (1989), 
the KS statistic is 0.122 and the x 2 statistic is 5.0. This 
indicates that the observed PNLF is unlikely to be drawn 
from the extrapolated one at the 90% confidence level. 

Improved simulations of the empirical PNLF were pre- 
sented for both the [Om]5007A and H/3 lines by Mendez 




Apparent Hp 4861A mag (corrected for foreground reddening) 



Figure 16. The H/3 planetary nebulae luminosity function for the 
central 25 deg 2 of the LMC, derived from a homogeneous sample 
of LMC PNe with rri486i <27.0. The data have been binned into 
0.2 mag intervals and plotted in log space. Error bars are based 
on poisson statistics. The solid line is the truncated exponential 
curve as predicted by Ciardullo et al. (1989) convolved for errors 
and placed at the best fit position. It marks out magnitude 16.8 
as the position of m* . The curve matches the brightest 6 magni- 
tudes of the function quite well. It is assumed that the falloff at 
magnitudes rri486i > 24 is due to incompleteness. 



and Soffner (1997). By varying the estimated fraction of 
stellar ionising luminosity absorbed by the nebula, better 
known as the absorbing factor fj,, and using available evolu- 
tionary tracks to produce a good representation of varying 
T e ff, mass and log L for 3000 100 year intervals between 
and 30,000 years, they produced good representations to the 
standard predictions of Ciardullo (1995). Their Fig. 6 simu- 
lation of the H/3 PNLF shows a peak in the distribution at 
precisely the same position of M = 3.8 found in this study 
using real LMC PNe. This remarkable match to model pre- 
dictions proves that not all PNe on cooling tracks can have 
fl = 1, since that would create a strong hump at the peak 
and make fainter PNe extremely hard to find. The values of 
/x most likely vary between 0.1 and 1 between ages of 1 and 
30,000 years. 



5.1 PNe with low central star temperatures 

It is important to note that we have found 3 PNe (SMP 
numbers 64, 31, 55) with very high H/3 fluxes belonging to 
excitation class 1, the lowest category (Reid & Parker, 2009). 
In RPb, SMP64 and SMP31 were classified as possible PNe 
belonging to a class known as Very Low Excitation (VLE) 
PNe. They have central stars with T e jj so low (31,200 K to 
38,500 K) that the [O iii]5007A is extremely low by compari- 
son to H/3 (details in RPb). A further 4 PNe (SMP67, 77, 76 
and MG62) also have [O m]5007/H,3 flux ratios ~3, low ex- 
citation classes from 2 to 3 and T e // between 46,500 K and 
64,200 K. The presence of VLE PNe has previously been re- 
ported by Mendez et al. (1993) and can be seen in Fig. 4a 
of Doptia et al. (1992) for several Magellanic Cloud PNe. 

In order for the H/3 luminosity to remain constantly 
high in these PNe, the number of ionising photons from the 
central star must increase by a factor of 2.5 as the tempera- 
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Apparent Hp 4861 A mag (corrected for internal extinction) 

Figure 17. The bright end of the H/J planetary nebulae lumi- 
nosity function to mag 21.4 for the central 25 deg 2 of the LMC, 
binned into 0.2 mag intervals and plotted in linear space. The 
broken line is the truncated exponential curve as predicted by 
Ciardullo et al. (1989). Error bars are based on poisson statistics. 
The curve matches the bright end of the function well. This pro- 
vides the first, direct estimation of H/3 M* for the LMC using the 
bright end of the LMC PNLF. 



ture increases from T e // 30,000 K to 70,000 K. If these neb- 
ulae were optically thick, we would expect a proportional 
increase in the H/3 luminosity with increasing central star 
temperature. To keep these nebulae constantly bright on 
the horizontal portion of the post-AGB evolutionary tracks, 
Mendez et al. (1993) suggest they must be optically thin 
in the H Lyman continuum. The fraction of stellar ionising 
luminosity absorbed by the nebula, defined as the absorb- 
ing factor, n, (Mendez et al. 1992) may also play a defining 
role here. Another possibility is that many of these central 
stars are burning He rather than H. This might account for 
a decrease in luminosity as T e ff increases. 



5.2 [Oin]5007 vs U/3 magnitudes 

Comparisons have been previously drawn between the 
[Om]5007A and Ha + [Nil] line ratios in a given popu- 
lation of extragalactic PNe (eg. Ciardullo et al. 2002, 2004). 
The Ha, being the brightest line in the Balmer series, is 
the easiest hydrogen line to measure in extra-galactic PNe. 
Unfortunately, due to filter responses, it usually includes 
the adjacent [Nll]6548A + 6583A lines. The LMC is close 
enough that the H/3 line is easy to measure using our spec- 
troscopic observations. Since it is cleanly separated and yet 
relatively close to the [O in] 5007 line in the optical spectrum, 
the ratios will be less affected by errors in flux calibration 
and extinction estimates. 

The shape of the H/3 PNLF has strong similarities to 
that of [O ill] Figure [6] but there some differences. Once 
again there is a peak to the distribution, however, it is not 
as strongly marked as its counterpart in [O ill] . With the 
H/3 peak occurring at magnitude 22.4 (absolute mag 2.6), 
5 magnitudes below the brightest, 48% or 287 PNe have 
fluxes occurring at the bright end. This may be compared 
with the [O ill] function where the peak occurs at magnitude 
20.4 (absolute mag 0.4), however only 242 or 41% of PNe 



are found at the bright end of the peak. This is despite the 
[O ill] peak occurring 5.2 magnitudes below the brightest, 
which is relatively 0.2 magnitudes fainter than the peak in 
H/3. In general terms, across a broad range of evolution, this 
suggests that the H/3 flux increases and fades at a steadier 
rate than [O ill] flux. The [O ill] flux appears to be highly 
sensitive to evolutionary spikes and dips whereas the H/3 flux 
is more sensitive to the opacity of the nebula. It also sug- 
gests that the central star evolves over the dynamic lifetime 
of a PN as it ionises the uniformly expanding shells (Henize 
& Westerlund, 1963), affecting temperature, luminosity and 
excitation of the nebula. These parameters will be closely 
examined in a further paper in this series. 

All the previously known and new PN in this survey are 
included in Figure [18] where the [O ill] 5007A magnitudes 
are plotted against the H,9 magnitudes, both de-reddened 
and corrected for extinction. The correction between the 
line fluxes is generally very small. Even for a huge red- 
dening of E(B-V)=1.0, the differential reddening between 
[Om] 5007A and H/3 is less than 0.06 dex (Schild, 1977). 

Figure [18] shows that the [O ill] 5007A line intensities 
span ~10 magnitudes. The H/3 line intensities on the other 
hand only span ~8.5 magnitudes. A faint dotted line in- 
cluded in Figure[T8]represents the H/3 function, moved to the 
best fitting position with the [Om]5007 function. Statisti- 
cally, on a bin-to-bin basis, there is an agreement of 82% be- 
tween them. This, however, should be viewed only in terms 
of overall shape of the PNLF since [O m]5007/H/3 in individ- 
ual PNe may vary by as much as 2.5 magnitudes. With this 
in place, it becomes apparent that [O ill] 5007 fluxes extend 
approximately 1 magnitude brighter than the bright end of 
the H/? function. Assuming an exponential distribution, the 
KS statistic is 0.10 and the x 2 statistic is 9.72 indicating a 
92% probability that the two functions do not correspond. 
Nonetheless, we cannot statistically exclude the possibility 
that the two samples are essentially the same shape. 

At the bright end, the H/3 PNLF like that of 
[O ill] 5007A starts to climb almost immediately without any 
strongly defined dips but a series of peaks and troughs. This 
may indicate that the central star's progression onto the WD 
cooling track has little immediate impact on stellar luminos- 
ity. Although the stellar temperature begins to fall, causing 
a dip in the [Om]5007A line, the H/3 flux is maintained 
longer, probably due to the higher ionisation potential of H. 

Although the transition probability (3 to 1) between the 
[Om]5007A and 4959A lines (Osterbrock & Ferland, 2006) 
remains constant, our measurements of PNe in the LMC 
show that the H/3 flux varies widely by comparison to these 
lines. Combined [O iii]5007A + 4959A fluxes lie between 62 x 
and 0.9 x the strength of H/3 PN fluxes in the LMC survey 
sample. Figure [191 directly compares the measured fluxes. A 
fitted trend line shows a steady regression towards the line 
of equality with fainter magnitudes. PNe at the brightest 4 
magnitudes have consistently high 5007/4861 A ratios. The 
marked lack of low 5007/4861A ratios in this region suggests 
that these PNe are young, optically thick, medium to high 
excitation PNe with central stars which have not yet reached 
the WD cooling track. 

PNe with low [Om] fluxes, -0.9 to 2 times H/3 are 
found in both faint, highly evolved PNe and within mag- 
nitude ranges fainter than 18. Other diagnostics such as 
high [Nll]6583/Ha together with morphology, when avail- 
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Figure 18. The planetary nebulae luminosity functions of H/3 
and [O III] 5007, corrected for extinction, are directly compared 
using all PNe uncovered in the central 25 deg 2 of the LMC. The 
data have been binned into 0.2 mag intervals and plotted ac- 
cording to their measured and corrected magnitudes. The bright 
cut-off for H/3 occurs at magnitude 16.8, 2.4 magnitudes below 
that of [O in] 5007. At the peak of the distributions (shown) the 
Amagnitude (H/3- [O III] 5007) is reduced to 1.8 magnitudes and 
remains at approximately this level to the faint cut-off. 
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Apparent magnitude [O III] 5007 

Figure 19. A direct comparison of [O iii]5007A and H/3486lA ap- 
parent magnitudes corrected for internal extinction shows a 
steady regression towards the line of equality with fainter mag- 
nitudes. The faint trend line marks out this regression while the 
darker line represents the line of equality where the 5007A and 
486lA magnitudes would be equal. PNe in the brightest 4 mag- 
nitudes mainly exhibit high 5007/486lA ratios. 



able from HST images, help to distinguish these candidates 
as possible Very Low Excitation (VLE) PNe. The mean ratio 
([O m]5007,4959/rL9) across the entire central 25deg 2 survey 
region is 8.93±7.31. By far, the majority (94%) of [O m]5007 
+ 4959 fluxes are below 20 x H/3 fluxes. 

The H/3 bright cutoff occurs 2.4 magnitudes fainter 
than the equivalent position for [Om]5007A. Figure [19] 
shows that PNe with the brightest H/3 fluxes also have 



strong [Oin]5007A fluxes. In fact, sampling the brightest 
50 H/3 fluxes, the average [Om]5007A flux is 8.9 times H/3. 
These include mostly medium and some high excitation PNe 
with central star temperatures between T e // =31,800 K and 
136,000 K. 

The H/3 PNLF shows the same rather flat distribu- 
tion seen across the initial 2 brightest magnitudes of the 
[Om]5007A PNLF. The initial peak in the H/3 distribution 
is only 1.8 magnitudes fainter than the equivalent position 
for [Oin]5007A. Similarly the faintest PNe in each wave- 
length are generally 1.8 magnitudes in separation. Both of 
these PNLF estimates agree closely with the average magni- 
tude difference of 1.7 ±0.9 by directly comparing 591 LMC 
PNe. By comparison, if we average the magnitude difference 
between [O m]5007A and H/3 for the brightest 30 PNe in H/3, 
the result is a much higher 2.3 ±0.6 mag difference. The flat 
distribution at the bright end of the [Om]5007A PNLF is 
therefore the result of a bright extension, peculiar to the 
bright end of the [Om]5007A PNLF. This extension is an 
average 0.6 mag greater than the average difference for the 
remaining 95% of PNe. These results, despite the evolu- 
tionary effects, indicate that the temperature of the central 
star is the most dominant contributor to the PN luminosity 
function (eg. Jacoby 1989), giving rise to both the H/3 and 
[Om]5007A luminosity. 

The raw cumulative luminosity functions derived from 
the [O in] and H/3 magnitudes are shown in Figure 1201 In 
each case, the bright end of the functions are fitted with 
a straight line, indicating the upper luminosity cutoff. This 
line traditionally represents the observational basis for dis- 
tance estimates once corrections for reddening and luminos- 
ity efficiency of the [O ill] or H/3 line photons are made. 
The graphs clearly show however, that most PNe follow the 
steeper slope which represents the greatest growth, peak and 
fall in the function. The number of PNe on this line is equal 
in [O m]5007A and H/3 (70-75% for both). The stronger peak 
in the [O iii]5007A function results in the oscillation seen ei- 
ther side of the arrow, to the right and left of the central 
line. 

Tables 2] and [S] in the appendix provide a compendium 
of and [Om]5007A raw flux estimates for previously 
known and newly discovered RP PNe respectively in the 
Large Magellanic Cloud. The tables also include the mea- 
sured values of cH/3 and the estimated errors. Table [5] in- 
cludes the PN object probability where 'T' represents 'True', 
'L' represents 'Likely' and 'P' represents 'Possible'. 

It is anticipated that once the RP survey is extended 
to the outer LMC, a similar shape for the PNLF to that 
newly obtained for the central 25 deg 2 would be found. The 
increase in total numbers should create a much tighter fit 
to the exponential curve, providing a very accurate distance 
estimate. The combination of previously known and newly 
discovered PNe in this study currently provide the most ac- 
curate estimate for the shape of the LMC PNLF obtained 
thus far. 



6 CONCLUSION 

Due to the depth of the UKST multi-stacked Ha survey, the 
subsequently uncovered large RP PN sample is able to show 
the influence of population effects on the shape of the PNLF 
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Figure 20. The cumulative PNLF using the [O III] 5007 line (left) and H/34861 line (right). The bright end of the PNLF, traditionally 
used for distance determination, is marked on both wavelengths. Both lines equally contribute to this section of the function, making 
either a good choice for modeling the bright end and distance estimations. The stronger peak in the [Om]5007A function results in the 
oscillation seen either side of the arrow, to the right and left of the central line. 



for the first time. We present details of a new technique for 
the flux calibration of fibre-based spectra. Using this tech- 
nique, we provide carefully calibrated [Oiii]5007A and H/3 
flux estimates from 2dF and FLAMES spectra together with 
fluxes from calibrated longslit spectra for 584 LMC PNe. 
These fluxes are converted to magnitudes and used to con- 
struct a PNLF in [O in] and in H/3. The PNe are separated 
into previously known and newly discovered PNe where a 
significant increase in the number of LMC PNe below an 
apparent magnitude of 18.6 in [Oiii]5007 is possible thanks 
to this work. The [Oiii]5007A PNLF is used to directly es- 
timate the distance to the LMC. By fitting an exponential 
curve to the bright end of the function over a 6 magnitude 
range, an m* value of 14.02 is found. Applying the accepted 
M* value of -4.44, we find a distance modulus of 18.46 for 
the LMC. This is in excellent agreement with the best recent 
calibrations (see table 2). The bright end cutoff is found to 
be robust to small samples of bright PNe since our signifi- 
cantly increased PN samples have no effect on this fiducial. 

We find a peak in the [O ill] 5007 PNLF ~6 mags below 
the brightest. This peak was previously predicted but not 
seen in a real sample until now. The steep rise to the peak 
in the function over a 2 mag range was also not seen before 
but was also predicted due to the position of the WD cooling 
track and a mixture of H and He-burning stars in the LMC. 
Not withstanding incompleteness, the faint end of the PNLF 
can now be examined for the first time. Over a 4 mag range 
we find a linear drop-off in the number of PNe per mag bin. 

A comparison of the foreground de-reddened [O ill] 5007 
PNLF with the same [O ill] 5007 PNLF after correction for 
internal extinction reveals almost no change in the overall 
shape of the PNLF. If the observed PNLF is shifted by ap- 
plying 1 magnitude of extinction to each object, then the 
de-reddened PNLF can be almost exactly recovered. This 
remarkable finding means that although the PN extinction 
ranges from < c(H/3) < 1.3, from a global prospective, the 
PNLF acts as if each one has c(H/3) = 0.4. This will clearly 
make it easier to model the PNLF in distant galaxies, as it 
implies that, for a 1 dimensional analysis, there is only a 
mean offset to be applied. 

We have constructed the first LMC PNLF in H/3 using 



591 measured PN spectra. We find a value of 16.8 for the 
m 486i bright cut-off. This leads to an estimate of M* = -1.68 
for H/3-based PNLF distance estimates. 

A clear relation is found between the shape of the 
[O ill] 5007A and the H/3 PNLF with a mean difference of 
2.8 magnitudes at the bright cut-off decreasing to 1.8 mag- 
nitudes by the peak of the distributions or 6 magnitudes 
below mJoQ7' 

The new LMC PNLF is compared to a new PNLF con- 
structed using local Galactic PNe to 3kpc and the newest 
sample from the SMC. The magnitude range is almost ex- 
actly the same for both the LMC and local surveys. A small 
dip ~3 and ~4 mag below the brightest is evident in the lo- 
cal and LMC surveys respectively. The rollover at the faint 
end of the PNLF is generally interpreted as a sign of in- 
completeness in the survey, where the rollover occurs across 
the faintest 1 or 2 magnitudes. Using 584 PNe in the RP 
sample from this survey, the rollover occurs over 1.5 magni- 
tudes, causing the shape of the faint end to inversely reflect 
the shape of the PNLF at the bright end. This would indi- 
cate that the peak distribution found in this work is close to 
the true peak distribution for LMC PNe. This occurs at an 
absolute magnitude of ~0.8. By comparison, the peak dis- 
tribution of local Galactic PNe occurs at an absolute mag- 
nitude of ~2, 1.2 magnitudes fainter than that found for the 
LMC. 

Once we complete our survey by including the outer 
LMC, we expect the estimated extra 80 PNe at the bright 
end of the function will fill out much of the exponential 
curve. It should thereby clarify the position of the rise, 4 
mags below the brightest and the position of any dip if they 
exist. The extra objects are not expected to alter the dis- 
tance modulus of 18.46; similar to the result found by Ja- 
coby (1989). Final results will largely depend on individual 
luminosities up to 6 magnitudes below m*. 
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Table 4. A compendium of H/3 and [Oiii]5007A raw flux esti- 
mates for Large Magellanic Cloud PNe. * = A flux measured di- 
rectly from long-slit spectra, f = PNe with extremely high Ha/H/3 
>10 ratios which may have been affected by power-law photoion- 
isation or shock- heating (see section [3.31 1. 
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RP 


Common 


cH/3 


Flux H/3 


Flux [Oiii] 


Ref. 


Nom. 




4868A 

(org cm- 2 s- 1 ) 


5007A 

(erg cm" 2 


1800 


MG8 


0.38 


1.52E-14±7.6E-16 


1.24E-13 ±6.2E-15 


1801 


MG10 


0.43 


1.25E-14 ±6.2E-16 


1.34E-13 ±6.7E-15 


1802 


MG15 


0.34 


4.46E-15 ±4.4E-16 


1.78E-14±1.5E-14 


1894 


SMP27 


0.26 


1.72E-14 ±8.6E-16 


2.24E-13 ±1.1E-14 * 


1895 


MG9 


0.39 


9.07E-15 ±7.7E-15 


7.32E-14 ±3.6E-15 


1896 


MG11 


0.41 


1.23E-14 ±6.1E-16 


5.88E-14 ±3.2E-14 


1897 


MG12 


0.41 


9.59E-15 ±8. IE- 15 


4.42E-14 ±2.4E-14 


1898 


Mol7 


0.28 


5.89E-15 ±5.0E-15 


3.97E-14 ±3.3E-14 



Table 5. A compendium of K/3 and [Oiii]5007A raw flux esti- 
mates for RP PNe in the Large Magellanic Cloud. P = Probabil- 
ity where T = True, L = Likely, P = Possible. A detailed review 
of T, L and P is given in [RPb]. * = A flux measured directly from 
long-slit spectra, f = PNe with extremely high Ha/H/3 >10 ra- 
tios which may have been affected by power-law photoionisation 
or shock- heating (see section [3.3l l. 



RP P cH/3 Flux HyS Flux [Oiii] 

Rcf. 4861A 5007A 

(erg cm- 2 S- 1 ) (erg cm" 2 a ~ 1 ) 



Explanation of abbreviations used: J: Jacoby (1980), LM2: Lindsay & Mul- 
lan (1963), MG: Morgan & Good (1992), Mo: Morgan (1994), N: Henize 
(1956), SMP: Sanduleak ct al. (1978), Sa: Sanduleak (1984), WS: Wcst- 
erlund & Smith (1964) 
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8.47E-14 ±4.24E-15 


624 


T 





.41 


6 


.56E- 


16 ±6 


56E- 


17 


1.53E-15 ±1.30E-15 


980 


T 


0. 


38 


2 


.05E- 


13±1 


03E- 


14 


9.10E-13 ±4.55E-14 


637 


T 





.45 


2 


31E- 


15±1 


16E- 


16 


6.46E-15 ±3.23E-16 


992 


P 


0. 


43 


9 


.96E- 


14 ±4. 


98E- 


•15 


5.83E-14 ±2.91E-15 


645 


P 





.69 


1 


56E- 


13 ±7 


.82E- 


15 


6.18E-14±3.09E-15 


1006 


T 


0. 


35 


3 


.24E- 


15 ±1. 


62E- 


16 


2.84E-14 ±1.42E-15 


666 


T 





.41 


1 


.41E- 


14±7 


■ 03E- 


16 


6.33E-14±3.16E-15 


1012 


T 


0. 


80 


1 


57E- 


14 ±7. 


86E- 


•16 


5.25E-14 ±2.62E-15 


671 


T 





70 


1 


■ 08E- 


15 ±5 


41E- 


17 


1.24E-14±6.18E-16 


1018 


P 


0. 


57 


2 


.48E- 


15 ±1. 


24E- 


16 


4.46E-15 ±2.23E-16 


676 


T 





20 


1 


47E- 


15 ±7 


35E- 


17 


3.50E-15±1.75E-16 


1024 


T 


0. 


47 


1 


.82E- 


15 ±9 


09E- 


■17 


4.10E-15 ±2.05E-16 


678 


T 





.28 


2 


.87E- 


15±1 


.43E- 


16 


1.22E-14±6.11E-16 


1037 


T 


0. 


47 


4 


.44E- 


15 ±2. 


22E- 


•16 


1.60E-14 ±8.01E-16 
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HP P cH/3 Flux H/3 Flux [O III] RP P cH/3 Flux H/3 Flux [O III] 

Ref. 4868A 5007A Ref. 4868A 5007A 

(erg cm -2 s~ 1 ) (erg cm~ 2 s~ 1 ) (erg cm -2 s — 1 ) (erg cm -2 s~ 



1038 


T 


0. 


19 


2 


30E-15 ±1 


15E- 


•16 


6.72E-15 ±3.36E-16 


1040 


T 


0. 


38 


1 


97E-15 ±9 


86E- 


•17 


1.52E-14 ±7.61E-16 


1045 


L 


0. 


04 


2 


96E-16 ±2 


96E- 


•17 


1.45E-15 ±1.23E-15 * 


1046 


T 


0. 


35 


3 


53E-16 ±3 


53E- 


•17 


1.45E-15 ±1.24E-15 


1065 


P 


0. 


49 


5 


79E-15 ±2 


89E- 


•16 


1.67E-14 ±8.37E-16 * 


1066 


L 


0. 


28 


4 


69E-15 ±2 


35E- 


16 


1.16E-14 ±5.80E-16 * 


1067 


T 


0. 


43 


1 


30E-15 ±6 


48E- 


•17 


4.27E-15 ±2.14E-16 


1068 


T 


0. 


67 


6 


98E-16±6 


98E- 


•17 


2.73E-15 ±1.37E-16 


1069 


P 


0. 


42 


3 


39E-14 ±1 


69E- 


•15 


6.90E-14 ±3.45E-15 * 


1071 


P 






2 


58E-13±1 


29E- 


14 


8.94E-15 ±4.47E-16 * 


1072 


T 


0. 


26 


3 


26E-15 ±1 


63E- 


•16 


2.23E-14 ±1.12E-15 


1078 


T 


0. 


36 


3 


06E-15 ±1 


53E- 


•16 


4.11E-14 ±2.06E-15 


1080 


T 


0. 


44 


2 


.75E-15 ±1. 


37E- 


•16 


1.14E-14 ±5.70E-16 


1081 


T 


0. 


39 


2 


80E-15 ±1 


40E- 


•16 


1.94E-14 ±9.72E-16 


1084 


T 


0. 


56 


3 


51E-16 ±3 


51E- 


•17 


1.45E-15 ±1.23E-15 


1088 


T 


0. 


03 


4 


44E-17 ±6 


67E- 


18 


4.89E-16 ±4.89E-17 


1089 


T 


1. 


21 


3 


90E-15 ±1 


95E- 


•16 


3.60E-14 ±1.80E-15 * 


1090 


T 


0. 


48 


5 


92E-16 ±5 


92E- 


•17 


7.84E-15 ±3.92E-16 


1092 


T 


0. 


34 


9 


76E-15 ±4 


88E- 


•16 


1.99E-14 ±9.93E-16 


1093 


L 


0. 


45 


1 


17E-15 ±5. 


85E- 


•17 


1.84E-15 ±9.19E-17 


1095 


T 


0. 


44 


3 


41E-16 ±3 


41E- 


•17 


3.38E-15 ±1.69E-16 


1106 


L 


0. 


55 


1 


17E-14 ±5. 


83E- 


•16 


9.36E-14 ±4.68E-15 


1148 


T 


0. 


23 


3 


.88E-15 ±7 


63E- 


17 


4.00E-15 ±1.30E-15 


1168 


T 


0. 


31 


7 


22E-15 ±3 


61E- 


•16 


4.40E-14 ±2.20E-15 


1179 


T 


0. 


19 


1 


59E-15 ±7 


94E- 


•17 


4.07E-15 ±2.03E-16 


1183 


T 


0. 


42 


3 


62E-15 ±1 


81E- 


16 


2.73E-14 ±1.37E-15 


1184 


T 


0. 


40 


5 


86E-16 ±5 


86E- 


•17 


2.85E-15 ±1.43E-16 


1185 


L 


1. 


07 


1 


26E-15 ±6 


30E- 


•17 


1.17E-14 ±5.84E-16 * 


1186 


T 


0. 


00 


1 


27E-15 ±6 


36E- 


17 


1.19E-14 ±5.93E-16 


1188 


T 


0. 


32 


7 


70E-15 ±3 


85E- 


•16 


4.68E-14 ±2.34E-15 


1189 


T 


0. 


38 


3 


00E-16 ±1 


50E- 


• 16 


9.72E-15 ±4.86E-16 


1190 


T 


0. 


36 


4 


54E-15 ±2. 


27E- 


■16 


3.25E-14 ±1.63E-15 


1191 


T 


0. 


22 


3 


19E-14 ±7 


64E- 


16 


1.81E-13 ±7.44E-15 


1196 


T 


0. 


06 


1 


19E-14 ±2 


11E- 


•16 


4.15E-14 ±1.08E-15 


1197 


T 


0. 


23 


1 


39E-14 ±1 


06E- 


•16 


8.43E-14 ±1.04E-15 


1201 


L 


1. 


56 


2 


81E-16±8 


48E- 


17 


2.95E-15 ±3.17E-16 * 


1207 


T 


0. 


74 


5 


77E-16±5 


77E- 


•17 


7.76E-15 ±3.88E-16 


1208 


T 


0. 


44 


9 


76E-15 ±1 


62E- 


•16 


2.45E-14 ±5.42E-16 


1233 


L 


0. 


40 


9 


28E-16 ±7 


89E- 


16 


1.14E-14 ±5.72E-16 


1240 


T 


0. 


32 


2 


79E-15 ±1 


39E- 


•16 


1.14E-14 ±5.72E-16 


1241 


L 


0. 


15 


2 


07E-14 ±4 


35E- 


•16 


2.78E-14 ±6.39E-16 


1242 


P 






7 


11E-17 ±1 


07E- 


17 


3.26E-16 ±3.26E-17 


1245 


P 






7 


76E-16±7 


76E- 


17 


4.04E-15 ±2.02E-16 


1246 


T 


0. 


55 


6 


23E-15 ±3 


11E- 


•16 


2.99E-14 ±1.50E-15 


1259 


T 


0. 


05 


2 


67E-15 ±1 


33E- 


•16 


5.89E-15 ±2.94E-16 


1267 


T 


0. 


28 


1 


90E-15 ±9 


50E- 


17 


2.19E-14 ±1.09E-15 


1270 


T 


0. 


02 


3 


92E-15 ±1 


96E- 


■16 


2.12E-14 ±1.06E-15 


1281 


P 


0. 


56 


9 


46E-16 ±8 


04E- 


•16 


1.33E-15 ±1.13E-15 


1284 


T 


0. 


27 


2 


91E-15 ±1 


45E- 


16 


1.72E-14 ±8.62E-16 


1288 


T 


0. 


51 


3 


14E-15 ±1 


57E- 


•16 


2.53E-14 ±1.27E-15 


1289 


T 


0. 


36 


6 


20E-15 ±3 


10E- 


•16 


4.78E-14 ±2.39E-15 


1296 


T 


0. 


10 


3 


24E-15 ±1. 


62E- 


•16 


5.90E-15 ±2.95E-16 


1300 


P 


0. 


88 


1 


22E-14 ±6 


09E- 


•16 


1.26E-13 ±6.32E-15 


1303 


T 


0. 


23 


1 


83E-15 ±9 


13E- 


•17 


1.23E-14 ±6.13E-16 


1304 


T 


0. 


33 


3 


05E-1S ±1 


52E- 


•16 


3.71E-14 ±1.85E-15 


1308 


T 


0. 


45 


3 


01E-15 ±1 


50E- 


•16 


9.06E-15 ±4.53E-16 


1309 


T 


0. 


29 


1 


50E-15 ±7 


50E- 


•17 


9.82E-15 ±4.91E-16 


1310 


T 


0. 


62 


8 


65E-16 ±8 


65E- 


•17 


5.44E-15 ±2.72E-16 


1314 


T 


0. 


28 


9 


99E-16±8 


49E- 


16 


5.69E-15 ±2.84E-16 


1315 


T 


0. 


54 


1 


03E-15 ±5 


14E- 


•17 


5.76E-15 ±2.88E-16 


1317 


T 


0. 


72 


3 


41E-16 ±3 


41E- 


•17 


1.77E-15 ±8.86E-17 


1323 


T 


0. 


38 


1 


02E-15 ±5 


10E- 


•17 


2.02E-15 ±1.01E-16 


1324 


T 


0. 


48 


2 


12E-1S ±1. 


06E- 


•16 


1.45E-14 ±7.25E-16 


1336 


T 


0. 


51 


4 


31E-1S ±2 


15E- 


■16 


1.31E-14 ±6.55E-16 


1337 


T 


0. 


54 


6 


86E-15 ±3 


43E- 


•16 


2.08E-14 ±1.04E-15 


1338 


T 


0. 


32 


1 


50E-14 ±7 


49E- 


16 


4.24E-14 ±2.12E-15 


1341 


P 


0. 


40 


1 


27E-15 ±6 


35E- 


•17 


3.00E-15 ±1.50E-16 


1345 


T 


0. 


40 


9 


79E-15 ±4 


89E- 


•16 


2.75E-14 ±1.37E-15 


1352 


T 


0. 


57 


1 


96E-14±9 


82E- 


16 


7.06E-14±3.53E-15 


1353 


T 


0. 


26 


2 


18E-15 ±1 


09E- 


•16 


7.48E-15 ±3.74E-16 


1354 


T 


0. 


58 


7 


24E-16 ±7 


24E- 


•17 


3.02E-15 ±1.51E-16 


1357 


T 


0. 


73 


1 


66E-15 ±8 


29E- 


17 


1.81E-14 ±9.07E-16 


1358 


T 


0. 


43 


5 


09E-15 ±2 


54E- 


•16 


4.25E-14 ±2.12E-15 


1371 


T 


0. 


58 


1 


10E-14 ±5 


50E- 


•16 


4.66E-14 ±2.33E-15 


1375 


P 


0. 


42 


1 


14E-14 ±5. 


68E- 


•16 


2.11E-14 ±1.05E-15 


1376 


T 


0. 


47 


4 


93E-1S ±2 


47E- 


16 


2.93E-14 ±1.46E-15 


1378 


T 


0. 


70 


4 


84E-1S ±2 


42E- 


16 


1.30E-14 ±6.52E-16 


1387 


T 


0. 


41 


7 


93E-15 ±3 


97E- 


•16 


6.31E-14 ±3.16E-15 


1397 


P 


0. 


38 


1 


41E-13 ±7 


07E- 


15 


1.82E-13 ±9.08E-15 



1415 


T 


0.47 


1.44E-15 ±7 


20E- 


17 


8.34E-15 ±4.17E-16 


1416 


T 


0.33 


2.24E-15 ±1 


12E- 


16 


1.56E-14 ±7.78E-16 


1418 


T 


0.26 


1.21E-15 ±6 


07E- 


•17 


5.74E-15 ±2.87E-16 


1426 


T 


0.36 


4.21E-15 ±2 


10E- 


16 


1.60E-14 ±8.01E-16 


1429 


T 


0.26 


3.38E-15 ±1 


69E- 


•16 


2.65E-14 ±1.33E-15 


1431 


L 


0.46 


1.08E-15 ±5 


42E- 


17 


3.19E-15 ±1.60E-16 


1432 


T 


0.30 


9.18E-16±7 


81E- 


•16 


2.24E-15 ±1.12E-16 


1438 


T 


0.66 


3.56E-15 ±1 


78E- 


■16 


3.60E-14 ±1.80E-15 


1440 


T 


0.57 


1.28E-15 ±6 


40E- 


•17 


1.48E-14 ±7.42E-16 


1443 


T 


0.39 


5.36E-15 ±2 


68E- 


•16 


1.13E-14 ±5.63E-16 ' 


1444 


T 


0.43 


6.36E-15 ±3 


18E- 


16 


4.37E-14 ±2 


18E-15 


1446 


T 


0.70 


2.55E-15 ±1 


28E- 


•16 


2.77E-14 ±1. 


39E-15 


1447 


T 


0.01 


1.58E-15 ±7 


90E- 


•17 


5.86E-15 ±2 


93E-16 


1456 


T 


0.28 


3.99E-15 ±1 


99E- 


•16 


2.44E-14 ±1. 


22E-15 


1462 


T 


0.44 


2.26E-15 ±1 


13E- 


•16 


9.24E-15 ±4 


62E-16 


1463 


T 


0.23 


3.88E-15 ±1 


94E- 


16 


1.27E-14 ±6 


33E-16 


1465 


T 


0.80t 


8.08E-16 ±8 


08E- 


•17 


4.46E-15 ±2 


23E-16 


1466 


P 


0.50 


2.76E-14 ±1 


38E- 


•15 


3.46E-14 ±1 


73E-15 


1467 


T 


0.74 


6.31E-16±6 


31E- 


•17 


2.17E-15 ±1 


09E-16 


1474 


T 


0.01 


6.30E-16±6 


30E- 


•17 


5.52E-15 ±2 


76E-16 


1488 


T 


0.39 


7.67E-16±7 


67E- 


•17 


5.34E-15 ±2 


67E-16 


1502 


T 


0.37 


1.81E-14 ±9 


03E- 


•16 


9.16E-15 ±1 


58E-16 


1508 


L 


0.10 


2.95E-15 ±1 


47E- 


16 


2.32E-14 ±1 


16E-15 


1519 


T 


0.71 


1.60E-16 ±1 


60E- 


17 


1.31E-1S ±1 


11E-15 


1523 


T 


0.44 


3.89E-15 ±1 


94E- 


•16 


2.14E-14 ±1. 


07E-15 


1528 


T 


0.66 


8.42E-15 ±4 


21E- 


16 


8.23E-14 ±4 


12E-15 


1532 


T 


0.48 


1.25E-14 ±6 


27E- 


16 


8.01E-14 ±4 


00E-15 


1550 


T 


0.26 


2.23E-14 ±1 


11E- 


15 


8.59E-14 ±4 


29E-15 


1559 


P 


0.32 


1.11E-15 ±5 


57E- 


■17 


2.10E-1S ±1 


05E-16 


1579 


T 


0.37 


2.39E-15 ±1 


20E- 


•16 


1.14E-14 ±5 


69E-16 


1580 


T 


0.88 


3.49E-15 ±1 


75E- 


16 


1.04E-14 ±5 


22E-16 


1584 


T 


0.86t 


1.63E-15 ±8 


16E- 


•17 


8.21E-15 ±4 


11E-16 


1587 


T 


0.34 


7.41E-15 ±3 


70E- 


16 


7.02E-14 ±3 


5 IE- 15 


1595 


T 


0.12 


6.71E-15 ±3 


36E- 


■16 


8.20E-15 ±4 


10E-16 


1601 


P 


0.86 


8.78E-15 ±4. 


39E- 


■16 


1.18E-14 ±5 


92E-16 


1615 


P 


0.27 


8.66E-15 ±4 


33E- 


16 


1.31E-14 ±6 


57E-16 


1624 


T 


0.28 


5.98E-15 ±8 


53E- 


17 


1.53E-14 ±2 


91E-16 


1631 


T 


0.06 


3.37E-15 ±1 


69E- 


■16 


7.62E-15 ±3 


81E-16 


1634 


P 


0.91 


4.96E-15 ±6 


67E- 


■17 


1.50E-14 ±2 


84E-16 


1636 


T 


0.33 


5.62E-15 ±2 


81E- 


•16 


4.41E-14 ±2 


20E-15 


1638 


T 


0.19 


2.14E-15 ±1 


07E- 


•16 


1.12E-14 ±5.58E-16 ' 


1649 


P 


0.43 


6.95E-16 ±6 


95E- 


•17 


1.26E-15 ±1 


07E-15 


1659 


T 


0.31 


9.27E-16±7 


88E- 


16 


5.78E-15 ±2 


89E-16 


1660 


T 


0.30 


4.39E-15 ±2 


19E- 


16 


3.24E-14 ±1 


62E-15 


1664 


L 


0.46 


9.42E-16 ±8 


01E- 


•16 


2.49E-15 ±1 


25E-16 


1676 


T 


0.33 


1.24E-14 ±6 


20E- 


16 


5.08E-14 ±2 


54E-15 


1694 


T 


0.44 


1.74E-15 ±8 


72E- 


17 


1.44E-14 ±7. 


22E-16 


1695 


T 


0.36 


5.18E-15 ±2 


59E- 


16 


4.09E-14 ±2 


04E-15 


1696 


T 


0.41 


1.86E-15 ±9 


30E- 


17 


1.46E-14 ±7 


29E-16 


1697 


T 


0.50 


7.12E-16 ±7 


12E- 


•17 


3.57E-15 ±1 


78E-16 


1705 


T 


0.37 


1.50E-1S ±7 


48E- 


■17 


5.24E-15 ±2 


62E-16 


1709 


T 


0.40 


1.93E-15 ±9 


65E- 


•17 


1.92E-15 ±9 


62E-17 


1712 


P 


0.50 


1.04E-15 ±5 


21E- 


•17 


1.56E-15 ±1 


33E-15 


1714 


T 


0.46 


5.19E-15 ±2. 


60E- 


16 


3.25E-14 ±1 


63E-15 


1718 


T 


0.09 


5.87E-16 ±5 


87E- 


•17 


9.08E-16±9 


08E-17 


1719 


T 


0.52 


8.48E-16 ±8 


48E- 


•17 


1.44E-14 ±7. 


20E-16 


1720 


T 


0.48 


2.95E-15 ±1 


47E- 


•16 


2.07E-14 ±1 


04E-15 


1721 


T 


0.56 


5.11E-16 ±5 


11E- 


•17 


3.47E-15 ±1 


73E-16 


1725 


T 


0.31 


3.03E-15 ±1 


51E- 


•16 


2.05E-14 ±1 


03E-15 


1726 


T 


0.29 


1.34E-15 ±6 


69E- 


•17 


1.35E-14±6 


74E-16 


1727 


T 


0.23 


5.73E-15 ±2 


86E- 


16 


1.91E-14 ±9 


53E-16 


1731 


T 


0.30 


1.85E-15 ±9 


25E- 


•17 


8.28E-15 ±4 


14E-16 


1739 


T 


0.47 


1.30E-15 ±6 


50E- 


•17 


6.85E-15 ±3 


42E-16 


1740 


T 


0.27 


5.66E-15 ±2 


83E- 


16 


1.80E-14 ±9 


00E-16 


1742 


T 


0.21 


6.06E-16±6 


06E- 


•17 


1.95E-15 ±9 


77E-17 


1743 


T 


0.36 


4.58E-15 ±2 


29E- 


16 


2.99E-14 ±1 


49E-15 


1748 


T 


0.30 


4.04E-15 ±2 


02E- 


16 


2.70E-14 ±1 


35E-15 


1749 


T 


0.17 


9.26E-16±7 


87E- 


16 


7.58E-16±7. 


58E-17 


1753 


T 


0.33 


8.82E-16±8 


82E- 


17 


2.04E-14 ±1 


02E-15 


1756 


T 


0.37 


7.37E-16±7 


37E- 


•17 


3.08E-15 ±1 


54E-16 


1758 


T 


0.33 


5.66E-15 ±2 


83E- 


16 


2.39E-14 ±1 


20E-15 


1759 


T 


0.34 


4.31E-15 ±2 


16E- 


16 


1.45E-14 ±7. 


26E-16 


1771 


T 


0.58 


3.74E-15 ±1 


87E- 


•16 


1.26E-14 ±6 


31E-16 


1773 


T 


0.54 


5.88E-15 ±2 


94E- 


16 


1.46E-14 ±7 


31E-16 


1791 


T 


0.49 


1.06E-14±5 


29E- 


•16 


2.09E-14 ±1 


04E-15 


1796 


P 


0.84 


3.16E-15 ±1 


58E- 


16 


2.29E-15 ±1 


14E-16 


1803 


T 


0.31 


1.52E-15 ±7. 


62E- 


•17 


5.37E-15 ±2 


69E-16 


1805 


T 


0.43 


1.88E-15 ±9 


38E- 


•17 


7.51E-15 ±3 


76E-16 



28 Warren A. Reid and Quentin A. Parker 



HP 
Ref. 


P 


cH/3 


1807 


T 


0.21 


1808 


L 


0.11 


1811 


P 


0.99 


1813 


T 


0.19 


1819 


T 


0.60 


1822 


L 


0.43 


1823 


T 


0.42 


1835 


T 


0.27 


1848 


T 


0.41 


1853 


T 


0.51 


1862 


L 


0.62 


1864 


T 


0.31 


1868 


T 


0.38 


1872 


T 


0.42 


1876 


T 


0.25 


1878 


T 


0.51 


1886 


T 


0.07 


1900 


T 


0.28 


1904 


T 


0.11 


1921 


T 


0.48 


1922 


T 


0.42 


1930 


P 


0.09 


1934 


T 


0.31 


1938 


T 


0.18 


1954 


T 


1.17 


1955 


T 


0.46 


1956 


T 


0.78 


1957 


T 


0.32 


1958 


P 


0.73 


1962 


P 


0.06 


2193 


P 


0.34 



Flux H/3 
4868A 





(erg cm 2 


b- 1 


1 


1 


67E-15 ±8 


35E- 


•17 


3 


39E-15 ±1 


69E- 


•16 


8 


59E-16 ±8 


59E- 


•17 


9 


21E-16 ±7 


83E- 


•16 


2 


02E-16 ±1 


01E- 


•16 


9 


77E-1S ±4. 


88E- 


16 


8 


40E-16 ±8 


40E- 


•17 


5 


70E-15 ±2 


85E- 


•16 


1 


,57E-14±7. 


87E- 


•16 


3 


71E-15 ±1 


86E- 


•16 


7 


34E-16±7 


34E- 


•17 


1 


35E-14±6 


74E- 


•16 


4 


20E-15 ±2 


10E- 


•16 


1 


25E-14 ±6 


25E- 


•16 


2 


22E-1S ±1. 


11E- 


•16 


6 


07E-15 ±3 


04E- 


16 


1 


68E-15 ±8 


42E- 


•17 


8 


13E-15 ±4 


06E- 


•16 


5 


60E-15 ±2 


80E- 


•16 


4 


69E-15 ±2 


35E^ 


•16 


2 


38E-15 ±1 


19E- 


•16 


3 


26E-16±3 


26E- 


•17 


5 


20E-15 ±2 


60E- 


16 


1 


10E-14 ±5 


52E- 


•16 


8 


78E-16±8 


78E- 


•17 


5 


08E-15 ±2 


54E- 


16 


2 


14E-1S ±1. 


07E- 


•16 


4 


82E-16 ±2 


41E- 


•16 


3 


72E-16±3 


72E- 


17 


4 


96E-14 ±2 


48E- 


•15 


1 


32E-14 ±6 


59E- 


•16 



Flux [Om] 
5007A 

(org cm -2 s -1 ) 



4.68E-15 ±2.34E-16 
8.74E-15 ±4.37E-16 * 
1.59E-15 ±1.35E-15 
7.10E-15 ±3.55E-16 
1.86E-14 ±9.28E-16 
2.57E-14 ±1.28E-15 
5.68E-15 ±2.84E-16 
3.44E-14 ±1.72E-15 
4.59E-14 ±2.29E-15 
1.46E-14 ±7.32E-16 
3.52E-15 ±1.76E-16 
2.24E-14 ±1.12E-15 
5.91E-14 ±2.96E-15 
1.01E-13 ±5.03E-15 * 



4 


■ 20E-15 ±2 


10E- 


•16 


3 


01E-14 ±1 


50E- 


15 


1 


25E-14 ±6 


26E- 


16 


1 


80E-14 ±8 


98E- 


16 


5 


32E-14 ±2 


66E- 


15 


1 


,75E-14±8. 


74E- 


•16 


1 


92E-14 ±9 


60E- 


16 


1 


10E-15 ±9 


36E- 


16 


1 


53E-14 ±7 


66E- 


16 


5 


18E-14 ±2 


59E- 


15 


4 


82E-15 ±2 


41E- 


•16 


2 


04E-14 ±1 


02E- 


15 


4 


44E-14 ±2 


22E- 


15 


1 


34E-14 ±6 


68E- 


16 


1 


17E-15 ±9 


93E- 


16 


2 


00E-14 ±1 


00E- 


•15 


1 


.42E-14 ±7. 


11E- 


16 



